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ABSTRACT 

We present 33 GHz photometry of 103 galaxy nuclei and extranuclear star-forming complexes taken 
with the Green Bank Telescope (GBT) as part of the Star Formation in Radio Survey (SFRS). Among 
the sources without evidence for an AGN, and also having lower frequency radio data, we find a median 
thermal fraction at 33 GHz of «76% with a dispersion of ~24%. For all sources resolved on scales 
<0.5kpc, the thermal fraction is even larger, being >90%. This suggests that the rest-frame 33 GHz 
emission provides a sensitive measure of the ionizing photon rate from young star-forming regions, 
thus making it a robust star formation rate indicator. Taking the 33 GHz star formation rates as a 
reference, we investigate other empirical calibrations relying on different combinations of warm 24 yum 
dust, total infrared (IR; 8 — 1000 /im), Ha line, and far-UV continuum emission. The recipes derived 
here generally agree with others found in the literature, albeit with a large dispersion that most likely 
stems from a combination of effects. Comparing the 33 GHz to total IR flux ratios as a function of 
the radio spectral index, measured between 1.7 and 33 GHz, we find that the ratio increases as the 
radio spectral index flattens which does not appear to be a distance effect. Consequently, the ratio of 
non-thermal to total IR emission appears relatively constant, suggesting only moderate variations in 
the cosmic-ray electron injection spectrum and ratio of synchrotron to total cooling processes among 
star-forming complexes. Assuming that this trend solely arises from an increase in the thermal fraction 
sets a maximum on the scatter of the non-thermal spectral indices among the star-forming regions of 
o-qNt < 0.13. 

Subject headings: galaxies: nuclei - Hii regions - radio continuum: general - stars: formation 



1. INTRODUCTION 

^ ■ Radio emission from galaxies is powered by a combi- 
CD . nation of distinct physical processes. And although it is 
\^ • energetically weak with respect to a galaxy's bolometric 
■ luminosity, it provides critical information on the massive 
, star formation activity, as well as access to the relativis- 
tic [magnetic field -f cosmic rays (CRs)] component in 
the interstellar medium (ISM) of galaxies. 
] Stars more massive than ~ 8 Mq end their lives as 
core-collapse supernovae, whose remnants are thought 
to be the primary accelerators of CR electrons (e.g., 
iKovama et al.ll 995') giving rise to the diffuse synchrotron 
k> . emission observed from star-forming galaxies (jCondonl 
'jjj ■ [1991 . These same massive stars are also responsible for 
d [ the creation of Hii regions that produce radio free-free 
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emission, whose strength is directly proportional to the 
production rate of ionizing (Lyman continuum) photons. 

Microwave frequencies, nominally spanning '-^l — 
100 GHz and observable from the ground, are partic- 
ularly useful in probing such processes. The non- 
thermal emission component typically has a steep spec- 
trum {Si, oc where a ~ 0.8), while the thermal 
(free- free) component is relatively flat {a ^ 0.1). Ac- 
cordingly, for globally integrated measurements of star- 
forming galaxies, lower frequencies (e.g., 1.4 GHz) are 
generally dominated by non-thermal emission, while the 
observed thermal fraction of the emission increases with 
frequency, eventually being d ominated by free-free emis- 
sion once beyond -30 GHz (Co ndon fc Yinl [l99l . For 
typical Hii regions, the thermal frac tion at 33 GHz can b e 
considerably higher, being ~80% (|Murphv et al.|[2011f) . 
Thus, observations at such frequencies, which are largely 
unbiased by dust, provide an excellent diagnostic for the 
current star formation rate (SFR) of galaxies. 

It is worth noting that the presence of an "anoma- 
lous" emission component in excess of free-free emis- 
sion between —10 and 90 GHz, generally attributed to 
electric dipole rotational emission from ultrasmall (a < 
10~^ c m) grains (e.g., 'Erickson''1957': 'Dra ine k, LazarianI 
Il998air b: Planck Collaboration et al. 20lif or magnetic 
dipole emission from thermal fluctuatio ns in the mag- 
netiz ation of interstellar dust grains (Draine fc LazarianI 
Il999f ). may comphcate this picture. For a single outer- 
disk s tar-forming region in NGC6946, IMurphv et aP 
(|2010f) reported an excess of 33 GHz emission relative 
to what is expected given existing lower frequency radio 
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data. This result has been interpreted as the first de- 
tection of so-cahed "anomalous" dust emission outside 
of the Milky Way. Given that the excess was only de- 
tected for a single region, this emission component may 
be negligible for globally integrated measurements. 

Due to the faintness of galaxies at high (i.e., >15 GHz) 
microwave frequencies, existing work has been restricted 
to the brightest objects, and small sample sizes. For 
example, past studies demonstrating the link between 
high-frequency free-free emission and massive star for- 
mation incl ude investigations of Galact ic star-forming re- 
gions fe.g.. iMezger fc Hende rson'1967 ), nea rby dwarf ir- 
regular ga laxies (e.g., i Kleiirfc Graevall986l ). galaxy nu- 
clei ( e.g.. [Turn er & Ro 19831 1994 V nearb y starbursts 
(e.g., iKlein et al.iil988t iTurner fc Holll985[ ). and super 
star clusters within nearby blue compact dwarf s (e.g., 
ITurner et al1[l99l iKobulnickv fc JohnsonI [199^ . And 
while these studies focus on the free-free emission from 
galaxies, each was conducted at frequencies <30GHz. 
With recent improvements to the backends of exist- 
ing radio telescopes, such as the Caltech Continuum 
Backend (CCB) on the Green Bank Telescope (GBT) 
and the Wideband Interfcrometric Digital ARchitecture 
(WIDAR) correlator on the Karl G. Jansky Very Large 
Array (VLA) , the availability of increased bandwidth is 
making it possible to conduct investigations for large 
samples of objects at frequencies '-^30 GHz. 

Here we present such an investigation using GBT 
33 GHz photometry for 103 galaxy nuclei and extranu- 
clear star-forming regions in a sample of 46 nearby 
galaxies included in the Star Formation in Radio Sur- 
vey (SFRS; see H2.ip . These galaxies, which are in- 
cluded in t he Spitzer Infrared N earby Galaxies Sur- 
vey f SINGS: iKennicutt et el !'2003') and Key Insights on 
Nearby Gala xies: a Far-Infra red Survey with Herschel 
(KINGFISH: IKennicutt et al.f201l[ ) legacy programs, are 
well studied and have a wealth of ancillary data available. 
We are currently in the process of collecting complemen- 
tary interfcrometric observations of these same targets 
using the VLA, which will be presented in a forthcom- 
ing paper. This paper is organized as follows: In §2 we 
describe our sample selection and the data used in the 
present study. In §3 we describe our analysis procedures. 
Our results are presented in §4 and discussed in §5. Fi- 
nally, in §6, we summarize our main conclusions. 

2. SAMPLE AND DATA 

In this section we describe the sample selection. We 
additionally present the GBT observations and provide 
a description of the ancillary data utilized for the present 
study. 

2.1. Sample Selection 

The Star Formation in Radio Survey (SFRS) sam- 
ple comprises nuclear and extranuclear star-forming re- 
gions in 56 nearby galaxie s (d < 30 Mpc) observed as 
part of the SINGS (Kcnn icutt et all 120031) and KING- 
FISH (IKennicutt et all 120 ll " legacy programs. Each of 
these nuclear and extranuclear star-forming complexes 
have ^ 1' X Q'.5 sized mid- infrared (i.e., low resolution 
from 5 — 14 /j,m and high resolution from 10 — 37 /im) 
spectral mappings carried out by the IRS instrument on 
board Spitzer, and 47" x 47" sized Herschel/PACS far- 
infrared spectral mappings for a combination of the prin- 



cipal atomic ISM cooling lines of [01] 63 /im, [OIII]88 /xm, 
[Nil] 122,205 /im, and [CII]158/im. Two galaxies that 
are exceptions include NGC 5194 and NGC 2403; these 
galaxies were part of the SINGS sample, but are not 
formally included in KINGFISH. They were observed 
with Herschel as part of the Very Nearby Galaxy Survey 
(VNGS; PI: C. Wilson). Similarly, there are additional 
KINGFISH galaxies that were not part of SINGS, but 
have existing Spitzer data: NGC 5457 (MlOl), IC342, 
NGC 3077, and NGC 2146. 

SINGS and KINGFISH galaxies were chosen to cover 
the full range of integrated properties and ISM con- 
ditions found in the local Universe, spanning the full 
range in morphological types, a factor of ^ 10^ in in- 
frared (IR: 8 — 1000 fiui) luminosity, a factor of ~ 10^ 
in Lia/Lopt, and a large range in star formation rate (< 
10~^ — 10 Mq yr~^). Similarly, spectroscopically targeted 
extranuclear sources included in SINGS and KINGFISH 
were selected to cover the full range of physical condi- 
tions and spectral characteristics found in (bright) in- 
frared sources in nearby galaxies, requiring optical and 
infrared selections. Optically selected extranuclear re- 
gions were chosen to span a large range in physical prop- 
erties including extinction-corrected production rate of 
ionizing photons [QiH°) ~ lO'^^ - 10^^ g-i]^ metaUic- 
ity (~ 0.1 — 3Zq), visual extinction {Ay < 4mag), 
radiation field intensity (100-fold range), ionizing stel- 
lar temperature (Toff ~ 3.5 — 5.5 x 10'' K), and local 
H2/H1 ratios (< 0.1— > 10). A sub-sample of infrared- 
selected extranuclear targets were chosen to span a range 
in /^(8/im)//i,(24/im) and fna/ Ui8 fj,m) ratios. 

The total sample over the entire sky consists of 118 
star-forming complexes (56 nuclei and 62 extranuclear 
regions), 103 of which (46 nuclei and 57 extranuclear re- 
gions; see Tables [1] and [2l respectively) are observable 
with the GBT (i.e., having S > —30°). Galaxy mor- 
phologies, adopted distances, and optically-defined nu- 
clear types are given in Table [TJ In Figure [TJ we show 
the location of each target region on Spitzer 24 /im im- 
ages; the corresponding circle diameters are 25", which 
match the FWHM of our lowest resolution data (i.e., the 
beam of the GBT 33 GHz radio data) for the present 
multiwavelength study. 

2.2. GBT Observations and Data Reduction 

Observations in the Ka band (26 — 40 GHz) were taken 
using the Caltech Continuum Backend (CCB) on the 
100 m Robert C. Byrd Green Bank Telescope (GBT; 
iJewell fc Prestagj I2004D over a two year period span- 
ning 2009 March through 2011 January. The CCB si- 
multaneously measures the entire Ka bandwidth over 
four equally spaced frequency channels (i.e., 27.75, 31.25, 
34.75, and 38.25 GHz) and synchronously reads out and 
demodulates the beam-switched signal to remove atmo- 
spheric fluctuation and/or gain variations. The FWHM 
of the GBT beam across the full Ka band was typically 
«25"among our sets of observations (see H3.ll) . Given the 
range of distances to the sample galaxies, this projects 
to linear scales of 0.37 — 3.7 kpc. Reference beams used 
for sky subtraction are measured by nodding away 
from the source, and their positions are identified on top 
of Spitzer 24 /im images in Figure [T] 

Observations were made using an "On-the-Fly Nod" 
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Table 1 

Nuclear Source Positions and 33 GHz Photometry 







R.A. 


Docl. 




TvDe^ 


Dist.*^ 


Nuc. Type'^ 


*^33 GHz 




(J2000) 


(J2000) 




('Mnr'l 




(mjy) 


NGC 337 





59 


50.3 


- 7 


34 


44 


SBd 


19.3 


SF 


2.59 + 0.21 


NGC 628 


1 


36 


41.7 


4-15 


46 


59 


SAc 


7.2 




0.11 + 0.03 


NGC 855 


2 


14 


3.7 


4-27 


52 


37 


E 


9.73 


SF 


0.68 + 0.05 


NGC 925 


2 


27 


17.0 


+33 


34 


42 


SABd 


9.12 


SF 


0.50 + 0.04 


NGC 1266 


3 


16 


0.8 


— 2 


25 


37 


SBO 


30.6 


AGN 


8.45 + 0.45 


NGC 1377 


3 


36 


38.9 


—20 


54 


6 


SO 


24.6 




<0.44 


IC 342 


3 


46 


48.5 


-f68 


5 


45 


SABcd 


3.28 


SF(*) 


25.39 + 1.27 


NGC 1482 


3 


54 


39.5 


-20 


30 


6 


SAO 


22.6 


SF 


10.78 + 0.56 


NGC 2146 


6 


18 


37.7 


+78 


21 


24 


Sbab 


17.2 


SF(*) 


76.16 + 3.83 


NGC 2403 


7 


36 


50.0 


+65 


36 


3 


SABcd 


3.22 


SF(*) 


0.28 + 0.03 


Holmbergll 


8 


19 


13.3 


+70 


43 


8 


Im 


3.05 


0.69 + 0.04 


NGC 2798 


9 


17 


22.8 


+41 


59 


57 


SBa 


25.8 


SF/AGN 


7.82 + 0.52 


NGC 2841 


9 


22 


2.7 


+50 


58 


36 


SAb 


14.1 


AGN 


1.53 + 0.08 


NGC 2976 


9 


47 


15.3 


+67 


55 





SAc 


3.55 


SF 


0.49 + 0.07 


NGC 3049 


9 


54 


49.6 


+ 9 


16 


17 


SBab 


19.2 


SF 


1.56 + 0.13 


NGC 3077 


10 


3 


19.1 


+68 


44 


2 


lOpec 


3.83 


SF(*) 


7.55 + 0.40 


NGC 3190 


10 


18 


5.6 


+21 


49 


54 


SAap 


19.3 


AGN(*) 


0.93 + 0.08 


NGC 3184 


10 


18 


16.7 


+41 


25 


27 


SABcd 


11.7 


SF 


<0.18 


NGC 3198 


10 


19 


54.9 


+45 


32 


58 


SBc 


14.1 


SF 


0.72 + 0.11 


IC 2574 


10 


28 


48.4 


+68 


28 


1 


SABm 


3.79 


SF(*) 


0.62 + 0.05 


NGC 3265 


10 


31 


6.7 


+28 


47 


48 


E 


19.6 


SF 


1.05 + 0.14 


NGC 3351 


10 


43 


57.8 


+11 


42 


14 


SBb 


9.33 


SF 


4.86 + 0.27 


NGC 3521 


11 


5 


48.9 


- 


2 


6 


SABbc 


11.2 


SF/AGN(*) 


1.45 + 0.10 


NGC 3627 


11 


20 


15.0 


+12 


59 


30 


SABb 


9.38 


AGN 


2.70 + 0.16 


NGC 3773 


11 


38 


13.0 


+12 


6 


45 


SAO 


12.4 


SF 


0.95 + 0.07 


NGC 3938 


11 


52 


49.5 


+44 


7 


14 


SAc 


17.9 


SF(*) 


<0.19 


NGC 4254 


12 


18 


49.4 


+14 


24 


59 


SAc 


14.4 


SF/AGN 


1.60 + 0.12 


NGC 4321 


12 


22 


54.9 


+15 


49 


21 


SABbc 


14.3 


AGN 


4.34 + 0.24 


NGC 4536 


12 


34 


27.1 


+ 2 


11 


17 


SABbc 


14.5 


SF/AGN 


17.49 + 0.88 


NGC 4559 


12 


35 


57.7 


+27 


57 


35 


SABcd 


6.98 


SF 


0.70 + 0.05 


NGC 4569 


12 


36 


49.8 


+13 


9 


46 


SABab 


9.86 


AGN 


1.81 + 0.13 


NGC 4579 


12 


37 


43.6 


+11 


49 


1 


SABb 


16.4 


AGN 


23.38 + 1.17 


NGC 4594 


12 


39 


59.4 


-11 


37 


23 


SAa 


9.08 


AGN 


69.30 + 3.47 


NGC 4625 


12 


41 


52.4 


+41 


16 


23 


SABmp 


9.3 


SF 


0.36 + 0.02 


NGC 4631 


12 


42 


5.9 


+32 


32 


21 


SBd 


7.62 


SF(*) 


3.83 + 0.21 


NGC 4725 


12 


50 


26.6 


+25 


30 


6 


SABab 


11.9 


AGN 


0.29 + 0.03 


NGC 4736 


12 


50 


53.0 


+41 


7 


14 


SAab 


4.66 


AGN(*) 


2.86 + 0.16 


NGC 4826 


12 


56 


43.9 


+21 


40 


59 


SAab 


5.27 


AGN 


5.88 + 0.31 


NGC 5055 


13 


15 


49.2 


+42 


1 


49 


SAbc 


7.94 


AGN 


2.58 + 0.25 


NGC 5194 


13 


29 


52.7 


+47 


11 


43 


SABbcp 


7.62 


AGN 


6.25 + 0.39 


NGC 5457 


14 


3 


12.6 


+54 


20 


57 


SABcd 


6.7 


SF(*) 


<0.73 


NGC 5474 


14 


5 


1.3 


+53 


39 


43 


SAcd 


6.8 


SF(*) 


0.17 + 0.05 


NGC 5713 


14 


40 


11.3 


- 


17 


26 


SABbcp 


21.4 


SF 


6.34 + 0.35 


NGC 5866 


15 


6 


29.5 


+55 


45 


47 


SO 


15.3 


AGN 


6.21 + 0.32 


NGC 6946 


20 


34 


52.3 


+60 


9 


14 


SABcd 


6.8 


SF 


15.84 + 0.80 


NGC 7331 


22 


37 


4.1 


+34 


24 


56 


SAb 


14.5 


AGN 


1.42 + 0.11 



(mJy) 



<0.09 



0.15 + 0.03 



0.29 + 0.07 



0.16 + 0.05 



^ Morphological type taken from the NASA Extragalactic Database (NED; http:/ /ncdwww. ipac.caltech.edu). 
^ Rcdshift-indepciidciit distance taken from the list compiled bv lKennicutt et al.l <|2011D . except for the two 
non-KINGFISH galaxies NGC 5194 JCiardullo et alj[200l) and NGC 2403 I IFreedman et al.|[200]]) . 
^ Nuclear Type based optical spectroscopy: SF — Star-Forming; AGN — Non-thermal emission as given in 
Table 5 of lMoustakas et al.l II2010I) or (*) Table 4 of lHo et all U997D . 

^ Measured 33 GHz flux densities before applying a correction for the oversubtraction of emission due to 
reference nods falling on bright regions of the galaxy. Corrections were only applied to cases where the 
estimated loss is >15%. 



variant of a symmetric nodding procedure (i.e., double- 
differencing technique; Readhead et al. 1989) with data 
being cohected continuously through an entire observa- 
tion, including slews between beams. This procedure al- 
ternately places the source of interest in each of the two 
beams of the Ka-band receiver in an A/B/B/A pattern. 
Such a sequence is able to cancel means and gradients in 
the atmospheric or receiver noise with time. Each nod- 
ding cycle lasts «70s: an on-source dwell time of 10 s 
for each of the 4 phases, along with 10 s spent on the 
initial position acquisition, and 10 s for the 2 slews be- 
tween beams. Thus, «40 s is spent on source per nodding 
cycle. A detailed description on the performance of the 



CCB receiver, the data reduction pipeline, and error esti- 
mates is given in Mason et al. ( 2009) . In addition to any 
systematic errors, we assign a calibration error of 10% to 
the flux density measurement at each frequency channel. 

Our observing strategy was constructed to make the 
most efficient use of the telescope. Thus, given the large 
range in brightness among our targeted regions, we var- 
ied the time spent on source based on an estimate of 
the expected 33 GHz flux density using the Spitzer 24 /im 
maps. The number of nods spent on each source was cho- 
sen to reach a S/N ratio of «10, while never spending 
<5min (i.e., w4 nods) nor >70min (i.e., «60 nods) of 
total integration time on a given source. Naturally, not 
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Table 2 

Extranuclecir Source Positions and 33 GHz Photometry 



ID 






R.A. 


Dccl. 


oo vj rxz 


OO yjin^ 








(J2000) 


(J2000) 


(mjy) 


(mJy) 


NGC 628 Enuc. 1 


1 


36 


45.1 


+15 47 51 


0.54 + 0.17 




NGC 628 Enuc. 2 


1 


36 


37.5 


+15 45 11 


0.35 ± 0.04 




NGC 628 Enuc. 3 


1 


36 


38.8 


+15 44 25 


0.54 + 0.04 


0.46 + 0.04 


NGC 628 Enuc. 4 


1 


36 


35.5 


+15 50 11 


0.23 + 0.05 


<0.16 


NGC 2403 


Enuc. 


1 


7 


36 


45.5 


+65 37 


1.68 + 0.14 




NGC 2403 


Enuc 


2 


7 


36 


52.7 


+65 36 46 


1.33 + 0.13 




NGC 2403 


Enuc 


3 


7 


37 


6.9 


+65 36 38 


3.04 + 0.19 




NGC 2403 


Enuc. 


4 


7 


37 


17.9 


+65 33 46 


1.37 + 0.08 




NGC 2403 


Enuc. 


5 


7 


36 


19.5 


+65 37 4 


1.28 + 0.21 




NGC 2403 


Enuc. 


6 


7 


36 


28.5 


+65 33 50 


0.73 + 0.06 




NGC 2976 


Enuc. 


1 


9 


47 


7.8 


+67 55 51 


3.62 + 0.21 




NGC 2976 


Enuc. 


2 


9 


47 


24.1 


+67 53 56 


2.12 + 0.15 


1.92 + 0.15 


NGC 3521 


Enuc. 


1 


11 


5 


46.3 


- 4 9 


0.47 + 0.02 


0.25 + 0.02 


NGC 3521 


Enuc. 


2 


11 


5 


49.9 


- 3 38 


0.90 + 0.09 




NGC 3521 


Enuc. 


3 


11 


5 


47.6 


+ 33 


0.37 + 0.03 




NGC 3627 


Enuc. 


1 


11 




16.2 


+12 57 50 


2.42 + 0.15 




NGC 3627 


Enuc. 


2 


11 


20 


16.3 


+12 58 44 


6.16 + 0.32 




NGC 3627 


Enuc. 


3 


11 


20 


16.0 


+12 59 52 


0.77 + 0.08 


0.62 + 0.08 


NGC 3938 


Enuc. 


1 


11 


52 


46.4 


+44 7 


0.15 + 0.02 


0.10 + 0.02 


NGC 3938 


Enuc. 


2 


11 


53 


0.0 


+44 7 54 


<0.66 




NGC 4254 


Enuc. 


1 


12 


18 


49.1 


+14 23 58 


1.10 + 0.08 




NGC 4254 


Enuc. 


2 


12 


18 


44.6 


+14 24 24 


0.19 + 0.04 


<0.13 


NGC 4321 


Enuc. 


1 


12 


22 


58.9 


+15 49 35 


0.15 + 0.05 


<0.14 


NGC 4321 


Enuc. 


2 


12 


22 


49.8 


+15 50 29 


0.16 + 0.04 


<0.13 


NGC 4631 


Enuc. 


1 


12 


41 


40.8 


+32 31 50 


0.90 + 0.08 




NGC 4631 


Enuc. 


2 


12 


42 


21.3 


+32 33 6 


0.67 + 0.06 




NGC 4736 


Enuc. 


1 


12 


50 


56.2 


+41 7 19 


2.58 + 0.15 




NGC 5055 


Enuc. 


1 


13 


15 


58.0 


+42 25 


0.43 + 0.12 


<0.37 


NGC 5194 


Enuc. 


1 


13 


29 


53.1 


+47 12 39 


<0.87 




NGC 5194 


Enuc. 


2 


13 


29 


44.1 


+47 10 21 


1.72 + 0.24 




NGC 5194 


Enuc. 


3 


13 


29 


44.6 


+47 9 55 


1.00 + 0.23 




NGC 5194 


Enuc. 


4 


13 


29 


56.2 


+47 14 7 


2.06 + 0.23 


1.77 + 0.23 


NGC 5194 


Enuc. 


5 


13 


29 


59.6 


+47 14 


2.05 + 0.24 


1.82 + 0.24 


NGC 5194 


Enuc. 


6 


13 


29 


39.5 


+47 8 35 


0.92 + 0.07 


0.85 + 0.07 


NGC 5194 


Enuc. 


7 


13 


30 


2.5 


+47 9 51 


<0.64 




NGC 5194 


Enuc. 


8 


13 


30 


1.6 


+47 12 51 


0.97 + 0.27 




NGC 5194 


Enuc. 


9 






uy.y 


+47 11 12 


1.03 + 0.22 


0.78 + 0.22 


NGC 5194 


Enuc. 


10 


13 


29 


56.7 


+47 10 45 


0.85 + 0.22 


<0.66 


NGC 5194 


Enuc. 


11 


13 


29 


49.7 


+47 13 28 


0.69 + 0.10 


< 0.30 


NGC 5457 


Enuc. 


1 


14 


3 


10.2 


+54 20 57 


0.39 + 0.12 




NGC 5457 


Enuc. 


2 


14 


2 


55.0 


+54 22 26 


0.29 + 0.08 




NGC 5457 


Enuc. 


3 


14 


3 


41.3 


+54 19 4 


8.27 + 0.52 




NGC 5457 


Enuc. 


4 


14 


3 


53.1 


+54 22 5 


2.11 + 0.22 




NGC 5457 


Enuc. 


5 


14 


3 


1.1 


+54 14 27 


2.49 + 0.27 




NGC 5457 


Enuc. 


6 


14 


2 


28.1 


+54 16 26 


3.38 + 0.33 




NGC 5457 


Enuc. 


7 


14 


4 


29.3 


+54 23 45 


4.69 + 0.35 




NGC 5713 


Enuc. 


1 


14 


40 


12.1 


- 17 47 


1.18 + 0.14 




NGC 5713 


Enuc. 


2 


14 


40 


10.5 


- 17 47 


1.48 + 0.15 




NGC 6946 


Enuc. 


1 


20 


35 


16.7 


+60 11 


0.74 + 0.05 


0.39 + 0.05 


NGC 6946 


Enuc. 


2 


20 


35 


25.5 


+60 10 1 


2.36 + 0.16 




NGC 6946 


Enuc. 


3 


20 


34 


51.9 


+60 12 45 


1.01 + 0.07 




NGC 6946 


Enuc. 


4 


20 


34 


19.2 


+60 10 9 


3.01 + 0.15 


2.89 + 0.15 


NGC 6946 


Enuc. 


5 


20 


34 


39.3 


+60 4 54 


0.51 + 0.04 




NGC 6946 


Enuc. 


6 


20 


35 


6.1 


+60 11 


2.86 + 0.18 




NGC 6946 


Enuc. 


7 


20 


35 


11.2 


+60 9 


3.20 + 0.20 




NGC 6946 


Enuc. 


8 


20 


34 


32.5 


+60 10 22 


1.80 + 0.14 




NGC 6946 


Enuc. 


9 


20 


35 


12.7 


+60 8 53 


2.49 + 0.17 





^ Measured 33 GHz flux densities before applying a correction for the over- 
subtraction of emission due to reference nods falling on bright regions of the 
galaxy. Corrections were only applied to cases where the estimated loss is 
>15%. 
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all nods were deemed usable given weather conditions 
and technical difficulties, and sources were revisited if 
possible. We list the number of nods used and taken, 
along with the corresponding time spent on source, for 
the nuclear and extranuclear regions in Appendix Tables 
[6] and [71 respectively. 

2.3. Ancillary UV, Optical, Infrared, and Radio Data 

Galaxy Evolution Explorer (GALEX) far-UV (FUV; 
1528 A) and near-UV (NUV; 2271 A) data were taken 
from the GALEX archive (GR6) and will be included 
in the GALEX Large Galaxy Atlas (M. Seibert et al. 
2012, in preparation). We also made use of star masks 
generated by the GALEX team to identify and remove 
foreground stars for the analysis. The angular resolu- 
tions of the FUV and NUV images are 4'.'25 and 5'.'25, 
respectively, while the calibration uncertainty for these 
data is wl5% in both bands. 

The Hg imagin g used in the analysis is taken from 
iLerov et"aLl ()2012D , where details about the data quality 
and preparation (e.g., correction for [Nil] emission) can 
be found. Like the UV data. Ha images were corrected 
for foreground stars. The typical resolution of the Ha 
images is Ri2", and the calibration uncertainty among 
these maps is taken to be «20%. 

Archival Spitzer 24 /im data were largely taken from 
the SINGS and Local Volume Legacy (LVL) legacy pro- 
grams, and have a calibration uncertainty of ~5%. De- 
tails on the associated observati on strategies and data 
reduction steps c an be found in iDale et all ()2007D and 
iDale et all ()2009[ ). respectively. Two galaxies, IC 342 and 
NGC 2146, were not a part of SINGS or LVL; t heir 24 ^m 
imaging comes from lEngelbracht et al.l ()2008l ) . 

We additionally made use of Hersc hel 70, 100, 160, 
and 2 50 /im data from KINGFISH (see iKennicutt et al.l 
1201 If). Observation s with the PACS instrument 
(jPoglitsch et"all IMol) were carried out at 70, 100, 
and 160 /xm in the Scan-Map mode and reduced by 
the Scanamorp hos data reduction pipeline, version 12.5 
(jRoussell 120121 ). The goal of the Scanamorphos reduc- 
tion is to preserve low-surface-brightness emission. The 
250 iJj n observations were made using the SPIRE instru- 
ment (|GrifRn et al.ll2010D and reduced using the HIPE 
version spire-5.0.1894. For NGC 2403 and NGC 5194, we 
make use of Herschel 70, 160, and 250 /im imaging ob- 
tained as part of the VNGS program. Information on 
the obse rvational strategy and data processing can be 
found in lBendo et all (|201ll ) for NGC 2403 and E. Men- 
tuch (2012, in preparation) for NGC 5194. 

Ancillary radio data at 1.365 and 1.697 GHz for 24 of 
the sample galaxies are available from the Westerbork 
Synth esis Radio Telescope SINGS survey (|Braun et al.l 
l2007f ). The intrinsic FWHM of the radio beams is ap- 
proximately 11" east-west by 11/ sin ^" north-south at 
1.5 GHz and scales as wavelength, where 5 is the source 
declination. For 5 galaxies (i.e., NGC 628, NGC 3627, 
NGC 4254, NGC 4321, NGC 4569) the north-south axis 
of the WSRT-SINGS beam at 1.7 GHz is larger than 
median FWHM of the 33 GHz GBT beam (including 
the la scatter), ranging between 36'.'5 and 44". Simi- 
larly, for 7 galaxies (i.e., NGC 628, NGC 3627, NGC 4254, 
NGC 4321, NGC 4569, and NGC 4725, and NGC 4826) 
the north-south axis of the WSRT-SINGS beam at 



1.4 GHz is larger than the median -t-lcr scatter on the 
33 GHz GBT beam, ranging between 29" and 55'.'5. Since 
this will affect the accuracy of the matched photome- 
try with the GBT 33 GHz data, we do not include these 
sources in the radio spectral index analysis. The flux 
density calibration of the radio maps is better than «5%. 

3. PHOTOMETRY AND ANALYSIS 

In the following section we describe our procedure to 
produce resolution-matched photometry among all data 
sets. Special considerations when dealing with the GBT 
single-beam photometry are also described in detail. 

We also note that the beam of the GBT averages over 
significantly large physical areas given the distances of 
the sample galaxies (e.g., as large as ~3.7kpc with a 
median linear scale of Ri0.9kpc among all nuclei and ex- 
tranuclear sources targeted), and likely contains a num- 
ber of non-coeval Hii regions, photodissociation regions, 
and diffuse emission. This is similar to the physica l scale 
of wO.Skpc investigated bv iMurphv et al.l (|2011f ). who 
found that conducting their analysis with and without 
local background subtractions did not qualitatively af- 
fect their main conclusions. We therefore do not attempt 
to subtract local background estimates, which would add 
more uncertainty to the present analysis. 

3.1. GBT Observations: Additional Gonsiderations 

Among all sets of observations, we measure a median 
beam width over all channels of 25'.'07 with a dispersion 
of 3'.'1. However, given that we are targeting resolved 
sources, and the beam size varies over the full Ka band, 
we apply a correction to the flux densities at each fre- 
quency channel as if their beam we re at the nominal 25" . 
We use the scaling factors given in IMurnhv et all (fMol) . 
since the distance to NGC 6946 is close to the median dis- 
tance among the entire SFRS sample. These scale factors 
were derived by computing the photometry for NGC 6946 
on the 8.5 GHz map using the beam sizes for each fre- 
quency channel averaged over that entire observing ses- 
sion, and are 0.89, 0.95, 1.00, and 1.04 at 27.75, 31.25, 
34.75, and 38.25 GHz, respectively. While the exact val- 
ues of these corrections factors will change depending on 
source geometry/morphology, the final 33 GHz photome- 
try, averaged over the entire band, should be robust. The 
combination of resolved sources and the changing beam 
size across the Ka band will complicate the interpreta- 
tion of in-band radio spectral indices, and is therefore 
not presented. 

While the unblocked aperture and active surface of the 
GBT mitigate the significance of sidelobes, we inspected 
calibration scans from each observing session to deter- 
mine if there were instances of days exhibiting sidelobes 
with large amplitudes. The added power may cause an 
overestimation of the 33 GHz flux density. We find that 
sidelobes, when measurable, have an amplitude that is 
<2% of the beam peak, on average, over the entire 2 
year observing campaign. At this level, our results should 
not be significantly affecte d by t he presence of sidelobes 
given that Murph y et a ll (|2010l ) found negligible differ- 
ences in their ancillary radio, submillimeter, and infrared 
photometry using a model beam having sidelobes at 5% 
of the beam peak. 

The averaged Ka-band flux densities, weighted by the 
errors from each channel over the full band, are given 
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Table 3 

Nuclear Photometry from Ancillary Data 





f,,f24 um) 
^mjy; 


Frn /10~^1 

f OVfT C t f^TTl ^ ^ 
^Ci b Clll J 


•-^1 . 7GHz 
^mjy; 


c 

'-'1 .4GHz 
(mjy) 


/41528A) 
(mJy) 


/,(227lA) 
(mJy) 


(erg s ^ cm ^) 


E{B - y)^ 

(Uld.gj 


NGC 337 


241.0 ± 42.7 


23.6 ± 4.2 






1.7± 


0.4 


3.8± 


0.9 


8.3± 


2.2 


0.112 


NGC 628 


57.4 ± 10.2 


8.3 ± 1.6 


















0.070 


NGC 855 


60.8 ± 10.8 


7.6± 1.5 






1.0± 


0.2 


1.6 ± 


0.4 






0.071 


NGC 925 


46.2 ± 8.2 


7.3± 1.3 


1.8 ± 0.3 


2.1 ± 0.4 


2.2 ± 


0.5 


3.0 ± 


0.7 


3.6± 


0.9 


0.076 


NGC 1266 


811.6 ± 143.8 


75.2 ± 27.7 






0.02 ± 


0.00 


0.14± 


0.03 






0.098 


NGC 1377 


1804.0 ± 319.7 


64.0± 11.7 


















0.028 


IC 342 


9561.4 ±1694.3 


584.5 ±320.8 






19.6 ± 


4.4 


119.4 ±27.1 






0.559 


NGC 1482 


1116.7± 197.9 


152.8 ± 37.7 






0.09 ± 


0.02 


0.3± 


0.1 






0.040 


NGC 2146 


7414.2 ±1313.8 


892.0 ±164.1 


504.3 ±89.4 


584.1 ±103.5 


0.16± 


0.04 


0.5 ± 


0.1 






0.096 


NGC 2403 


59.9 ± 10.6 


11. 3± 2.3 


1.4± 0.2 


1.8± 0.3 


1.5± 


0.3 


2.5 ± 


0.6 


2.7± 


0.7 


0.040 


Holmberg II 


35.2 ± 6.2 


1.8± 0.3 


1.0 ± 0.2 


1.2 ± 0.2 














0.032 


NGC 2798 


956.1 ± 169.4 


120.7 ± 22.7 


















0.020 


NGC 2841 


44.9 ± 8.0 


4.4 ± 0.9 


3.7± 0.7 


4.4 ± 0.8 


0.6± 


0.1 


1.0± 


0.2 


1.0 ± 


0.3 


0.015 


NGC 2976 


93.6 ± 16.6 


10.1 ± 2.0 


1.9± 0.3 


2.1 ± 0.4 


1.2± 


0.3 


2.0± 


0.5 


4.4 ± 


1.2 


0.072 


NGC 3049 


387.3 ± 68.6 


20.6 ± 3.9 






1.3± 


0.3 


2.2 ± 


0.5 


3.9 ± 


1.0 


0.038 


NGC 3077 


868.6 ± 153.9 


72.5 ± 12.5 














21. 8± 


5.7 


0.067 


NGC 3190 


108.5 ± 19.2 


22.7 ± 4.0 






0.07 ± 


0.01 


0.3± 


0.1 


0.4 ± 


0.1 


0.025 


NGC 3184 


129.7 ± 23.0 


11.6 ± 2.0 


2.1 ± 0.4 


2.3 ± 0.4 


0.4± 


0.1 


0.9 ± 


0.2 


0.8 ± 


0.2 


0.017 


NGC 3198 


364.6 ± 64.6 


20.8 ± 3.6 


2.9 ± 0.5 


3.2 ± 0.6 


0.2 ± 


0.1 


0.4 ± 


0.1 


0.7± 


0.2 


0.012 


IC 2574 


41.2 ± 7.3 


3.0 ± 0.5 


0.7± 0.1 


0.9 ± 0.2 


1.0± 


0.2 


1.2 ± 


0.3 


4.3 ± 


1.1 


0.037 


NGC 3265 


278.8 ± 49.4 


18.4 ± 3.2 






0.4 ± 


0.1 


0.8 ± 


0.2 






0.024 


NGC 3351 


1435.0 ± 254.3 


103. 7± 23.1 






2.0 ± 


0.4 


5.3 ± 


1.2 


11. 8± 


3.1 


0.028 


NGC 3521 


287.3 ± 50.9 


37.9 ± 9.0 






0.5 ± 


0.1 


2.0 ± 


0.4 


5.7± 


1.5 


0.057 


NGC 3627 


633.2 ± 112.2 


90.4 ± 16.2 






0.5 ± 


0.1 


1.7± 


0.4 


4.9 ± 


1.3 


0.033 


NGC 3773 


128.0 ± 22.7 


9.4 ± 1.8 






3.2 ± 


0.7 


4.0 ± 


0.9 






0.027 


NGC 3938 


68.2 ± 12.1 


11. 3± 2.1 


1.3± 0.2 


1.5 ± 0.3 


0.7 ± 


0.2 


1.4 ± 


0.3 


1.2± 


0.3 


0.021 


NGC 4254 


436.0 ± 77.3 


54.8 ± 16.2 






1.1± 


0.2 


3.1 ± 


0.7 


6.4± 


1.7 


0.039 


NGC 4321 


723.2 ± 128.2 


76.6 ± 14.3 






2.4 ± 


0.5 


5.6 ± 


1.3 


6.3± 


1.7 


0.026 


NGC 4536 


2474.1 ± 438.4 


193.3 ± 33.2 






0.3 ± 


0.1 


0.7 ± 


0.2 


7.0 ± 


1.8 


0.018 


NGC 4559 


84.2 ± 14.9 


12.0 ± 2.0 


2.5± 0.4 


3.0 ± 0.5 


1.2 ± 


0.3 


1.9 ± 


0.4 


4.0 ± 


1.0 


0.018 


NGC 4569 


701.4 ± 124.3 


47.3 ± 9.5 






0.9 ± 


0.2 


3.3 ± 


0.8 


8.9± 


2.3 


0.047 


NGC 4579 


170.7 ± 30.2 


15. 5± 2.7 






0.3 ± 


0.1 


0.8 ± 


0.2 


5.3± 


1.4 


0.041 


NGC 4594 


110.9 ± 19.6 


7.9 ± 1.6 






0.6 ± 


0.1 


1.5 ± 


0.3 






0.051 


NGC 4625 


45. 4± 8.1 


6.4 ± 1.2 






1.1 ± 


0.2 


1.6 ± 


0.4 


1.7± 


0.5 


0.018 


NGC 4631 


566.5 ± 100.4 


79.2 ± 13.9 


39.9 ± 7.1 


45. 8± 8.1 


1.6± 


0.4 


2.4 ± 


0.6 


4.8 ± 


1.3 


0.017 


NGC 4725 


54.2 ± 9.6 


4.4± 0.8 


0.3 ± 0.0 




0.17± 


0.04 


0.4 ± 


0.1 


0.3 ± 


0.1 


0.012 


NGC 4736 


1094.8 ± 194.0 


150.4 ± 33.1 


23.8 ± 4.2 


24.8 ± 4.4 


1.7± 


0.4 


5.8± 


1.3 


4.9 ± 


1.3 


0.018 


NGC 4826 


1088.9 ± 193.0 


167.6 ± 30.5 


33.9 ± 6.0 




1.0± 


0.2 


2.7± 


0.6 


17.1 ± 


4.5 


0.041 


NGC 5055 


353.0 ± 62.6 


76.8 ± 13.6 


18.1 ± 3.2 


20. 6± 3.7 


0.6± 


0.1 


1.9± 


0.4 


5.8 ± 


1.5 


0.018 


NGC 5194 


608.9 ± 107.9 


95.2 ± 26.0 


52.2 ± 9.3 


58.5 ± 10.4 


2.2 ± 


0.5 


6.1 ± 


1.4 


11. 1± 


2.9 


0.035 


NGC 5457 


177.7 ± 31.5 


20.6 ± 3.7 






0.9± 


0.2 


1.9± 


0.4 


2.0 ± 


0.5 


0.009 


NGC 5474 


14.4 ± 2.6 


2.5 ± 0.5 






1.2 ± 


0.3 


1.5 ± 


0.3 


1.7± 


0.5 


0.010 


NGC 5713 


1137.6 ± 201.6 


85.6 ± 16.6 






1.7± 


0.4 


3.6± 


0.8 


5.6± 


1.5 


0.039 


NGC 5866 


95.5 ± 16.9 


33.8 ± 5.9 






0.12 ± 


0.03 


0.6± 


0.1 






0.013 


NGC 6946 


4355.3 ± 771.8 


301.0 ± 60.0 


83.8 ±14.8 


91. 9± 16.3 


0.4± 


0.1 


2.0± 


0.5 


18.1 ± 


4.8 


0.343 


NGC 7331 


301.9 ± 53.5 


49.0 ± 8.6 


17.2 ± 3.0 


20.1 ± 3.6 


0.3± 


0.1 


1.6± 


0.4 


2.3 ± 


0.6 


0.091 



^ Galactic extinction taken from (ISchlegel et al.ll99Sfl used to correct the Ha and GALEX FUV and NUV flux densities assuming Ay / E{B — V) — 3.1 
and the modeled extinction curves of IWeingartner Sz Drainel l l20Qin ; IDraine) ]2Q03I) . 

in Tables [T] and [2] along with uncertainties; the corre- 
sponding effective frequency is given for each source in 
Appendix Tables [6] and [71 and is w33 GHz for each po- 
sition (i.e., a median of 32.92 GHz with a dispersion of 
0.79 GHz, corresponding to a ~2% scatter in flux den- 
sities assuming a spectral index near 33 GHz of '^0.5). 
For sources not detected at the 3 a level, we list the 3 a 
upper limit. Flux densities of the individual channels, 
before applying any corrections^ are given with 1 a errors 
in the Appendix Tables [Hand [7] for the galaxy nuclei and 
extranuclear regions, respectively. We do not give up- 
per limits for the individual channels, but rather list the 
actual measured values and estimated errors. 

Since the reference beam throw is only lf3, real signal 
in our reference positions is a concern; having reference 
positions on the galaxy may result in an underestimation 
of the 33 GHz flux density (see FigureH]). For NGC 6946, 
we note that the positions presented here correct those 



originally given in iMurphv et al.l ()2010( ). however this 
slight correction does not impact their results. To quan- 
titatively assess the severity of having real signal in the 
reference beam for our measurement, we estimate over- 
subtractions using the 24 /im data by comparing the flux 
densities measured at the reference positions compared 
to the on-source position. To correct for these losses, we 
add back to the 33 GHz flux densities the average flux 
densities measured among the reference positions using 
modeled 33 G Hz maps by scaling t he 24/.tm data using 
Equation 2 of lMurphv et all (|2006bD . The minimum flux 
density of the reference beam positions for each target 
was then taken as the local sky and subtracted. 

These corrections are typically small, having a median 
value of w2% among all sources. We conservatively apply 
corrections to only those 22 sources estimated to be miss- 
ing >15% of the actual flux density based on the 24 ^m 
photometry alone. The median correction to the 33 GHz 
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Table 4 

Extranuclear Photometry from Ancillary Data 



Galaxy 


/.(24/im) 


Fir/10-" 5i.7GHz 


5l.4GHz /.(1528A) 


/,.(2271 A) 


/hc«/10-13 E{B - vr 




(mjy) 


(ergs cm ^) (mJy) 


(mJy) (mJy) 


(mJy) 


(ergs -"-cm ^) (mag) 



NGC 628 Enuc. 1 


NGC 628 Enuc. 2 


NGC 628 Enuc. 3 


NGC 628 Enuc. 4 


NGC 2403 


Enuc 


1 


NGC 2403 


Enuc 


2 


NGC 2403 


Enuc 


3 


NGC 2403 


Enuc 


4 


NGC 2403 


Enuc 


5 


NGC 2403 


Enuc 


6 


NGC 2976 


Enuc 


1 


NGC 2976 


Enuc 


2 


NGC 3521 


Enuc 


1 


NGC 3521 


Enuc 


2 


NGC 3521 


Enuc 


3 


NGC 3627 


Enuc 


1 


NGC 3627 


Enuc 


2 


NGC 3627 


Enuc 


3 


NGC 3938 


Enuc 


1 


NGC 3938 


Enuc 


2 


NGC 4254 


Enuc 


1 


NGC 4254 


Enuc 


2 


NGC 4321 


Enuc 


1 


NGC 4321 


Enuc 


2 


NGC 4631 


Enuc 


1 


NGC 4631 


Enuc 


2 


NGC 4736 


Enuc 


1 


NGC 5055 


Enuc 


1 


NGC 5194 


Enuc 


1 


NGC 5194 


Enuc 


2 


NGC 5194 


Enuc 


3 


NGC 5194 


Enuc 


4 


NGC 5194 


Enuc 


5 


NGC 5194 


Enuc 


6 


NGC 5194 


Enuc 


7 


NGC 5194 


Enuc 


8 


NGC 5194 


Enuc 


9 


NGC 5194 


Enuc 


10 


NGC 5194 


Enuc 


11 


NGC 5457 


Enuc 


1 


NGC 5457 


Enuc 


2 


NGC 5457 


Enuc 


3 


NGC 5457 


Enuc 


4 


NGC 5457 


Enuc 


5 


NGC 5457 


Enuc 


6 


NGC 5457 


Enuc 


7 


NGC 5713 


Enuc 


1 


NGC 5713 


Enuc 


2 


NGC 6946 


Enuc 


1 


NGC 6946 


Enuc 


2 


NGC 6946 


Enuc 


3 


NGC 6946 


Enuc 


4 


NGC 6946 


Enuc 


5 


NGC 6946 


Enuc 


6 


NGC 6946 


Enuc 


7 


NGC 6946 


Enuc 


8 


NGC 6946 


Enuc 


9 



167.5 ± 29.7 
87.9 ± 15.6 

86.2 ± 15.3 
25.8 ± 4.6 

247.5 ± 43.9 
189.3 ± 33.5 
808.9 ±143.3 
68.8 ± 12.2 
171. 3± 

27.7 ± 
411. 4± 
186.5 ± 

23.8 ± 
138.9 ± 

30.3 ± 
405.9 ± 



30.4 

4.9 

72.9 

33.1 

4.2 

24.6 

5.4 

71.9 



1389.2 ±246.2 



139.3 ± 

46.8 ± 
57.4 ± 
114.4± 

49.7 ± 

54.9 ± 

47.8 ± 

46.8 ± 

84.7 ± 

506.4 ± 

78.3 ± 
250.9 ± 
299.9 ± 
227.1 ± 

168.5 ± 
196.5 ± 

66.9 ± 
180.8 ± 

327.5 ± 

146.1 ± 
175.7 ± 

78.4 ± 

97.8 ± 
60.8 ± 

1279.4 ±226.7 

143.3 ± 25.4 

186.6 ± 
248.0 ± 
162. 7± 
97.8 ± 
248.6 ± 

129.4 ± 
211. 9± 
92.6 ± 

102.2 ± 
36.4 ± 

445.3 ± 
454.0 ± 

233.4 ± 
367.6 ± 



24.7 

8.3 

10.2 

20.3 

8.8 

9.7 

8.5 

8.3 

15.0 

89.7 

13.9 

44.5 

53.1 

40.2 

29.9 

34.8 

11.9 

32.0 

58.0 

25.9 

31.1 

13.9 

17.3 

10.8 



33.1 
43.9 
28.8 
17.3 
44.1 
22.9 
37.5 
16.4 
18.1 
1.5 

78.9 
80.4 
41.4 
65.1 



10.0 ± 
5.5 ± 
6.4± 
2.1± 

20.5 ± 

18.3 ± 
45.8 ± 

7.9 ± 
14.0 ± 

3.3± 
32.8 ± 
18.8 ± 

2.2 ± 

19.7 ± 
3.7± 

32.6 ± 
111.7± 

19.6 ± 
6.2 ± 
3.4± 

10.0 ± 
6.5 ± 
7.1± 
6.5 ± 
6.3± 
8.8± 

54.4 ± 
8.0± 

29.1 ± 
37.4 ± 

22.8 ± 
21. 1± 

23.2 ± 
8.1± 

17.0 ± 
31. 2± 

17.3 ± 

24.4 ± 
9.8± 

12.3 ± 
5.4± 

50.2 ± 

10.9 ± 
10.9 ± 
17.8 ± 
8.1 ± 

14.2 ± 

20.3 ± 
11. 5± 
16.3 ± 

10.1 ± 
10.3 ± 
4.7± 
41.8± 

43.5 ± 
24.1 ± 
36.8 ± 



1.8 
1.0 
1.1 
0.4 
3.8 
5.0 
10.5 
1.4 
2.5 
0.6 
7.4 
3.5 
0.4 
3.6 
0.7 
5.8 
19.2 
3.4 
1.1 
0.6 
2.4 
1.3 
1.3 
1.2 
1.1 
1.6 
9.8 

I. 4 
8.8 

II. 4 
6.2 
4.1 
4.5 
1.6 
3.3 
9.0 
3.2 
5.9 
2.3 
2.3 
0.9 
9.9 
1.9 
2.2 
3.3 
1.9 
2.7 
3.7 
2.0 
2.9 
1.8 
1.8 

0.8 
7.4 
7.6 
4.3 
6.7 



4 


3± 





8 


4 


7± 





8 


3 


7± 





6 


4 


0± 





7 


9 


6± 


1 


7 


9 


9± 


1 


8 


2 


2± 





4 


2 


5± 





4 


3 


6± 





6 


3 


9± 





7 


1 


5± 





3 


1 


8± 





3 


5 


8± 


1 





6 


0± 


1 


1 


3 


2± 





6 


3 


3± 
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1.4± 0.2 
0.9 ± 0.2 



4.4 ± 

6.0 ± 

14.0 ± 
2.4 ± 

12.1 ± 

7.4 ± 
7.6 ± 
9.6 ± 
9.2 ± 
2.7± 
4.7± 
10.7± 

5.5 ± 

9.1 ± 
6.1 ± 



0.8 
1.1 
2.5 
0.4 
2.2 
1.3 
1.3 
1.7 
1.6 
0.5 
0.8 
1.9 
1.0 
1.6 
1.1 



3.7 ± 
6.5 ± 

3.4 ± 
3.3 ± 
2.2 ± 
11.8± 

9.8 ± 

6.5 ± 
8.5 ± 



0.7 
1.2 
0.6 
0.6 
0.4 
2.1 
1.7 
1.1 
1.5 



1.7± 0.3 
0.9± 0.2 



4.8 ± 
6.8± 
14.2 ± 

2.9 ± 
14.0 ± 
8.4± 
8.4± 
11.0± 
10.5 ± 
3.1± 
5.3 ± 
11.8± 
6.0 ± 
10.4 ± 
7.3 ± 



0.8 
1.2 

2.5 
0.5 

2.5 
1.5 
1.5 

1.9 

1.9 
0.6 
0.9 

2.1 
1.1 

1.8 
1.3 



4.3 ± 

7.4 ± 
3.8 ± 
3.8 ± 
2.7± 
13.0 ± 
10.9 ± 

7.5 ± 
9.6± 



0.8 
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0.7 
0.5 
2.3 
1.9 
1.3 
1.7 



0.6 ± 





1 


0.9 ± 


0.2 


2.2± 


0.6 


0.070 


0.5 ± 





1 


0.7 ± 


0.2 


1.3± 


0.3 


0.070 


1.5± 





3 


2.0± 


0.4 


2.4± 


0.6 


0.070 


0.4± 





1 


0.5 ± 


0.1 


1.2± 


0.3 


0.072 


6.4± 


1 


5 


8.2 ± 


1.9 


13.3 ± 


3.5 
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3.7± 
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0.040 
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1.1± 





3 


2.5 ± 


0.6 


7.3 ± 


1.9 


0.033 


1.8± 
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3.4± 
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3.0 ± 


0.8 


0.033 


0.7± 
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1.0± 


0.2 


1.1± 


0.3 
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0.5 ± 
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0.6± 


0.1 
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0.021 


0.7± 
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1.2± 
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0.5 ± 
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0.2 


1.3± 
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0.4± 
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0.8± 


0.2 


1.0 ± 
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0.026 


0.7± 
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1.1± 
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1.1± 


0.3 


0.026 


2.0± 
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4.5 ± 


1.2 
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1.3± 





3 


1.6± 
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3.9 ± 


1.0 
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4.1± 
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1.7± 





4 


2.9± 
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5.6 ± 


1.5 


0.036 


2.1 ± 





5 


3.4± 


0.8 


4.8 ± 


1.3 


0.036 


1.1± 





2 


1.7± 


0.4 


0.3 ± 


0.1 


0.035 


2.1± 





5 


3.4± 


0.8 


0.4 ± 


0.1 


0.035 


0.3± 





1 


0.4± 


0.1 


1.9 ± 


0.5 


0.038 


1.2± 





3 


1.8± 


0.4 


5.0 ± 


1.3 


0.036 


2.3± 





5 


3.5 ± 


0.8 


4.8 ± 


1.2 


0.035 


1.3± 





3 


2.3± 


0.5 


2.2± 


0.6 


0.035 


1.0± 





2 


1.8± 


0.4 


2.5± 


0.7 


0.036 


0.7± 





2 


1.2± 


0.3 


1.3± 


0.4 


0.035 


0.6± 





1 


1.1± 


0.2 


1.3± 


0.3 


0.009 


0.8± 





2 


1.0± 


0.2 


2.3± 


0.6 


0.009 


2.9± 
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3.5 ± 
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19.6 ± 
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0.009 


o.o m 


n 
u 


Q 


4 9-1- 




o.u in 


2 1 


n noQ 
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2.9± 





7 


3.3± 


0.8 


10.0 ± 


2.6 


0.009 


3.9± 





9 


4.3 ± 


1.0 


10.5 ± 


2.8 


0.009 
0.010 


0.4± 





1 


1.0± 


0.2 


1.1± 


0.3 


0.039 


0.6± 





1 


1.2± 


0.3 


2.2± 


0.6 


0.039 


2.2 ± 





5 


4.7± 


1.1 


10.5 ± 


2.8 


0.342 


6.2 ± 


1 


4 


11.4± 


2.6 


31.4± 


8.2 


0.344 


2.4± 





6 


4.7± 


1.1 


13.7± 


3.6 


0.342 


1.9± 





4 


4.1 ± 


0.9 






0.343 


1.4± 





3 


3.0± 


0.7 


8.9 ± 


2.3 


0.338 


2.9± 





7 


6.6± 


1.5 


23.2 ± 


6.1 


0.342 


2.5 ± 





6 


5.9± 


1.3 


22.3 ± 


5.9 


0.342 
0.343 


1.7± 





4 


4.1 ± 


0.9 


14.2 ± 


3.7 


0.342 



^ Galactic extinction taken from jSchlegel et aUl99SD used to correct the Ha and GALEX FUV and NUV flux densities assuming Av / E(B ~V) — 3.1 
and the modeled extinction curves o f lWeingartner fc Dranig 1 1200 1I') :[D ram 



flux density for these sources is «24%. The uncorrected 
33 GHz flux densities are given in the final columns of 
Tables ID and m We note that for the case of NGC 6946, 
using 8.5 GHz-derived 33 GHz maps resulted in correc- 
tions that agreed to those bas ed on the 24 /xm-de rived 
values to within a few percent (jMurphy et al.l[2"010t ). 

3.2. Ancillary Data 

The photometry was carried out on the ancillary data 
sets after matching their resolution, cropping each image 



to a common field-of-view, and re-gridding to a common 
pixel scale. To accurately match the photometry from 
these images to the GET measurements, maps were con- 
volved to the resolution of the GET beam in the Ka 
ba nd (i.e., 2 5") fol lowin g the image registration method 
of lAniano et al.l (|2011[ ) . Using the combination of all 
infrared data, total infrared (IR; 8 — 1000 fim) emission 
a nd uncertainty maps were constructed using the models 
of iDraine fcTi ()2007[ ) as described in G. Aniano (2012, 
submitted) . 
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For each data set, the flux density at the position of 
each GET pointing is the surface brightness multiphed by 
the beam sohd angle. Uncertainties on the photometry 
are estimated as a combination of the calibration and 
beam-size uncertainties; we assign a 17% uncertainty to 
account for the dispersion in the beam areas among our 
sets of GET observations as given in H'S.U 

In the case of the UV and Ha photome try, we correct 
each r egion for Milky Way extinction using lSchlegel et all 
(|T998l) assuming Ay/Ej B — V) = 3.1 and th e mod- 
eled extinction curves of IWeingartner fc Drain"i (|2001[ ): 
IDraind (|2003). The multi- wavelength photometry is 
given in Tables [3] and |31 along with the corresponding 
1-CT errors, for the nuclei and extranuclear regions, re- 
spectively. 

3.3. Star Formation Rate Calibrations 

For the non-AGN sources in the sample, we can esti- 
mate star formation rates using the new GET 33 GHz 
d ata. In this section w e present the calibrations given 
in lMurphv et al.l (|2011f) . where more details about their 
derivations can be fo und. These calibr ations, which up- 
date those found in iKennicuttI (I1998D and have been 
adopted by IKennicutt fc Evansl (120121), were calculated 
using Starburst99 fcei therer et al.lll999() for a common 
IMF so that each diagnostic can be compared fairly. A 
summary of these relations can be found in the Appendix 
Table E 

We choose a Kroupa (jKroupal I2001D IMF, having a 
slope of —1.3 for stellar masses between 0.1 — 0.5 Mq and 
—2.3 for stellar masses ranging between 0.5 — 100 Mq. 
Assuming solar metallicity and a continuous, constant 
star formation rate over ^100 Myr, Starburst99 stellar 
population models yield the following relation between 
the star formation rate and production rate of ionizing 
photons, Q{H°): 



SFR 

Mqyt- 



^ 7.29 X 10 



-54 



(1) 



Since the ionizing flux comes from very massive stars 
with lifetimes <10 Myr, we note that the coefficient in 
Equation [T] is nearly independent of starburst age under 
the assumption of continuous star formation so long as it 
is >10 Myr. Accordingly, such measurements sample the 
current (i.e., ~10 Myr) star formation activity. The in- 
tegrated UV spectrum is also dominated by young stars, 
however, it is sensitive to a significantly longer timescale 
of recent (~10 - 100 My r; lKennicuttl[T99l iCalzetti et all 
120051: iSalim et 311120071 ) star formation activity. We con- 
volve the output Starburst99 spectrum with the GALEX 
FUV transmission curve to obtain the following conver- 
sion between star formation rate and FUV luminosity. 



/ SFRfuv \ 



= 4.42 X 10" 



L 



FUV 



erg s 



(2) 



Similarly, we can write such an expression for the 
GALEX NUV band such that, 



/SFR 



•NUV 



VMoyr- 



= 7.15 X 10 



-44 



-t^NUV 

erg s^i 



(3) 



The ionizing photon rate can of course be expressed as 
a (extinction corrected) H recombination line flux, such 



that for Case E recombination, and assuming an elec- 
tron temperature Te = 10^ K, the Ha recombination line 
strength is related to the star formation rate by 



SFRhq 
Mq yr-1 



= 5.37 X 10" 



L 



Hq 



erg s 



(4) 



The above equation indicates that the star formation rate 
is directly proportional to the Ha line luminosity, assum- 
ing that a constant fraction (i.e., ~45%) of the ionized 
H atoms will emit an Ha photon as they recombine, and 
that the extinction correction is accurate. However, if 
a significant fraction of ionizing photons are absorbed 
by dust, the above equation will underestimate the star 
formation rate, and the = sign should be replaced by >. 

Similarly, at high radio frequencies, where r <g; 1, the 
ionizing photon rate is directly proportional to the ther- 
mal spectral luminosity, Lj, , varying only weakly with 
electron temperature To (jRubinlfl968[ ). such that 



= 6.3 X 10^^ 



104 K 



-0.45 



VGHz 



erg s" 



Hz- 



(5) 



Ey combining Equations [T] and [5l one can derive a rela- 
tion between the star formation rate and thermal radio 



emission: 



SFR^ 



Mq yr- 



= 4.6 X 10" 



lO^K 



-0.45 



(gHz) 



0.1 



erg s ^ Hz 



(6) 



As with the H recombination line fiuxes, Q{H'^), and con- 
sequently the star formation rate, may in fact be under- 
estimated by the free-free emission if a significant frac- 
tion of ionizing photons are absorbed by dust; in this 
case the = sign in the above equation should be replaced 
by >. However, it is worth noting that unlike free- free 
emission, which arises directly from the ionized gas, op- 
tical/NIR H recombination line fluxes may also suffer 
extinction internal to the Hii region itself, resulting in 
an even larger deflcit. For example, the extinction inter- 
nal to Hii regions in the starbursting dwarf galaxy NGC 
5253 is measured to be quite large {Ay = 16 — 18 mag; 
[Turner et all 12001 ). 

At lower radio frequencies, which are typically domi- 
nated by non-thermal synchrotron emission, calibrations 
between the supernova rate, and thus the star forma- 
tion rate, have been developed. From the output of 
Starburst99, which assumed a supernova cut-off mass of 
8 Mq, we find that the total core-collapse supernova rate, 
iVsN, is related to the star formation rate by. 



SFR 



M, 



©yr" 



86.3 



yr- 



(7) 



Work comparing the non-thermal spectral luminosity 
with the supernova rate in the Galaxy has yielded an 
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empirical calibration such that 



-NT 



erg s 



Hz- 



= 1.3 X 




GHz. 

(8) 

spectral index 
(jTammannl llQSa ICondon &: Yiijri990l ). By combining 
Equations [7] and |S1 we can express the star formation rate 
as a function of the non-thermal radio emission where 



where a^"^ is the non-thermal radio 



/ sfrn 



NT 



6.64x10 



-29 



VGHz/ 



L 



NT 



erg s ^ Hz 



VMgyr- 

(9) 

Since it will take '^30 Myr for 8 Mq stars to go supernova, 
and the radiating lifetimes of CR electrons can be on the 
order of tens of Myr in normal galaxies, the non-thermal 
emission is sensitive to slightly longer timescales than 
free-free or H recombination line emission. 

The observed radio continuum emission comprises both 
free-free and synchrotron emission. Therefore, we can 
combine Equations [6] and |9] to construct a single expres- 
sion for the star formation rate from the total radio con- 
tinuum emission at a given frequency such that: 



/ SFR^ 

\MQyr- 



10 



-27 



2.18 



15.1 



VlO^K 

1 



0.45 



VGHz 



(gHz) 
erg s~^ Hz~ 



(10) 



This equation essentially weights the observed radio con- 
tinuum luminosity based on the expected thermal frac- 
tion at a given frequency. As pointed out above, the 
thermal and non-thermal emission timescales are mis- 
matched, with free-free emission being sensitive to mas- 
sive stars with ages < 10 Myr and the non-thermal emis- 
sion being most sensitive to the lowest mass (i.e., >8 Mq) 
supernova progenitors having lifetimes of <30 Myr. How- 
ever, Equation [10] should hold under the assumption of 
continuous star formation on timescales >30 Myr. 

In our analysis we have assumed an electron tempera- 
ture of Tc — 10'' K when calculating star formation rates 
using Equations |6] and [101 However, variations in the 
actual electron temperatures will inject scatter into our 
observed trends. For instance, assuming a value of Tc 
as low as 5000 K will result in star formation rates that 
are 37 and 21% larger than assuming Tc = 10** K us- 
ing equations [6| and [TOl respectively. Thus, we expect 
any scatter introduced by the assumption of a constant 
electron temperature to be smaller than these values. 

4. RESULTS 

In the following section we estimate the fraction of 
thermal emission at 33 GHz and compare a number of 
star formation rate diagnostics with those derived using 
the GBT 33 GHz data. The assumption is made that 
the 33 GHz data provide the most reliable estimate as 
they are highly sensitive to free- free emission arising from 
the ionized gas in Hii regi ons. While we adopt the cal- 
ibrations given above Ci.e.. iMurphv et al.|[20ll . we de- 
rive our own empirical relations and compare with others 
found in the literature. 




0.1 



1.0 

s. 



10.0 
33GHz (mJy) 



100.0 



Figure 2. Radio spectral indices measured between 1.7 and 
33 GHz plotted against 33 GHz fiux densities for all 53 galaxy nu- 
clei (circles) and extranuclear star-forming regions (crosses) having 
1.7 GHz measurements. Nuclei identified as AGN in Table [l] are 
plotted using filled symbols. Upper limits for the 33 GHz flux den- 
sities, and corresponding spectral index measurements, are shown 
as arrows. Sources for which the major axis of the 1.7 GHz data 
was larger than the FWHM of the GBT beam at 33 GHz are not 
included due to imprecise matched photometry. 
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Figure 3. Estimates for the thermal fraction at 33 GHz plot- 
ted against the projected diameter of the GBT beam for all (i.e., 
41) non-AGN sources with 33 GHz detections and corresponding 
1.7 GHz data. Nuclei and extranuclear regions are shown by cir- 
cles and crosses, respectively. Only a weak trend is found with 
distance, suggesting that the 33 GHz star formation rates are not 
significantly affected by averaging over larger physical areas, which 
should increase the amount of non-thermal emission at 33 GHz, for 
the more distant sources. The seven sources having 33 GHz thermal 
fractions <50% are identified, five of which required a correction for 
badly placed off-nod positions. These are the same seven sources 
having star formation rate s discrepant by more than a factor of 
>1.5 in Figure |4] (see ^4.1|l . 

4.1. Radio Spectral Indices and Estimates for the 
Thermal Fraction at 33 GHz 

In Figure [2] the radio spectral indices measured be- 
tween 1.7 and 33 GHz are shown for all 53 sources hav- 
ing data at 1.7 GHz with a resolution that is equal to, or 
higher than, the GBT data. AGN are identified among 
the galaxy nuclei, typically having radio spectral indices 
that are steeper than the rest of the sources. Using the 
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Figure 4. Top: The 33 GHz star formation rates calculated in 
two ways for non-AGN sources liaving lower frequency radio data, 
allowing for an estimate for the 33 GHz thermal fraction. Star 
formation rates calculated using the thermal emission at 33 GHz 
and Equation [6] are given along the ordinate while star formation 
rates calculate d us ing the total 33 GHz spectral luminosity along 
with Equation 1101 assuming a non-thermal radio spectral index 
of a^T ^ _Q gg 

are given along the abscissa. Also shown is a 
one-to-one (dotted) line. The seven sources having star formation 
rates discrepant by more than a factor of >1.5 are the same seven 
sources identified in Figure [3] (see i|4.1l l. Bottom: The ratio of the 
two 33 GHz star formation rate estimates, illustrating the effect of 
assuming different values for the non-thermal radio spectral index, 
ranging between 0.6—1.05 (left to right) in increments of 0.05. The 
open circles assume a^'^ = —0.85, as in the above comparison. 



spectral index information over this rather long lever- 
arm, we can estimate the thermal emission fraction at 
33 GHz if we assume a fixed non-thermal spectral in- 
dex. We assume a constant non-thermal radio spectral 
index of a^''" = 0.85 given that this was the average non- 
thermal spectral index found amo ng the 10 star- f ormin g 
regions studied in NGC6946 bv iMurphv et all (|2011[ ). 
Furthermore , this is very close to the value found by 
iNiklas et all ()1997l i.e., a^'^ = 0.83 with a scatter of 
aQ,NT = 0.13) for a sample of 74 nearby galaxies. We do 
not estimate the thermal fractions for galaxy nuclei con- 
taining AGN since there is likely to be more variation in 
their i ntrinsic non-thermal slope. Following iKlein et al.l 
(jl984| ). and assuming that the free- free emission is opti- 
cally thin, we estimate the thermal fraction such that 
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Figure 5. The 33 GHz star formation rates normalized by the 
projected area of the GBT beam (i.e, spanning ^ 0.15 - llkpc^) 
plotted against distance for all (i.e., 81) non-AGN detected sources 
in the sample. Nuclei and extranuclear regions are shown by cir- 
cles and crosses, respectively. Only a weak trend is found with 
distance, suggesting that the 33 GHz star formation rates are not 
significantly affected by averaging over larger physical areas, which 
should increase the amount of non-thermal emission at 33 GHz, for 
the more distant sources. Error bars are not shown since they are 
typically smaller than the plotting symbols. 
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Figure 6. The 24 /im flux densities plotted against corresponding 
(corrected) 33 GHz flux densities for the entire sample of 103 galaxy 
nuclei (circles) and extranuclear star-forming regions (crosses). Nu- 
clei identified as AGN in Table [T] are indicated by filled symbols. 
Upper limits for the 33 GHz flux densities are shown as arrows. 
Error bars are not shown since they are typically smaller than 
the plotting symbols. Over-plotted is an ordinary least squares 
flt to the detected extranuclear regions and non-AGN nuclei (solid 
line), along with a one-to-one line scaled by the median 24fim to 
33 GHz flux density ratio < /^(24 ^m)/533 gHz >= 142 (dotted 
line). A number of AGN are clearly radio loud at 33 GHz relative 
to the main correlation between the 24 fim and 33 GHz flux densi- 
ties. The galaxy furthest above the flt to the main trend, which is 
detected at 24 /im but not at 33 GHz, is NGC 1377. This ga laxy is 
ident i fled a s being possibly a nascent starburst (e.g., IRoussel et al.l 
l200a . l20M) . 

where, for our specific case, i^i — 33 GHz, 1/2 = 1.7 GHz, 
a is the observed radio spectral index between and 1^2 , 
and we fix the thermal radio spectral index to 0.1. We 
find a median 33 GHz thermal fraction of «76% among 
all sources, with a lower and upper quartile of 64 and 
87%, respectively (see Figure [S]). While this method 
has been shown to overestimate thermal fractions for 
star-forming regions in a high-resolution radio study of 
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Figure 7. The ratio of 24 fim to 33 GHz star formation rates 
plotted against the 33 GHz star formation rates for all (i.e., 81) 
non-AGN detected sources in the sample. Nuclei and extranuclear 
regions are shown by circles and crosses, respectively. The scat- 
ter about unity is S!0.25dex. Extranuclear region 4 in NGC6946 
is identified using a square, as the 33 GHz flux density of this 
sourc es is complicated by the presence of anomalous dust emis- 
sion I IMurphv et al.|[20TOl: fScaife et al.ll20ia) . 



M33 (jTabatabaei et aLll2007[) . those authors assumed a 
much steeper non-thermal index (i.e., a^"^ = 1.0) than 
adopted here. Non-thermal spectral indices have been 
found to be significantly flatter than 1.0 around giant Hii 
regions based on comparisons between multifrequency 
radio and de-reddened Ha observations in NGC 6946 
(F. Tabatabaei et al. 2012, in preparation). Assuming 
Q-NT — X would increase our thermal fraction estimates 
to < f^^'^^^ >« 85 with a dispersion of 15%. 

In Figure H] we plot star formation rates for all 
non-AGN sources having low frequency radio data, al- 
lowing us to estimate the 33 GHz thermal fraction, 
in two ways. We first compute 33 GHz star forma- 
tion rates using only the thermal fraction at 33 GHz 
and Equation [51 In the top panel of 21 these are 
plotted against 33 GHz star formation rates using the 
total 33 GHz spectral luminosity with Equation [TOl 
again assuming a non-thermal spectral index of a^"^ = 
0.85. There are 6 sources (i.e, NGC 3938 Enuc. 1, 
NGC 4631, NGC 4631 Enuc. 2, NGC 5194 Enuc. 8, 
NGC 5194 Enuc. 10, and NGC 5194 Enuc. 11) with star 
formation rate estimates that are discrepant by more 
than a factor of ~2. This discrepancy may arise from 
our assumption for the non-thermal spectral index being 
too flat; these are the only sources shown in Figure [21 for 
which a?'^7GHz > 0.70, and are also the only sources in 
Figure Ohaving thermal fractions <40%. For the two ob- 
servations towards NGC 4631, we might expect a steeper 
spectral index as the galaxy is observed to be edge-on, 
and has a very prominent non-thermal radio component. 
Observing through the disk of this source could lead 
to a significant amount of diffuse non-thermal emission 
in the beam, steepening the observed spectrum. How- 
ever, for the remaining four sources, NGC 3938 Enuc. 1, 
NGC 5194 Enuc. 10, and NGC 5194 Enuc. 11, it is worth 
noting that the 33 GHz photometry required a correction 
for badly placed off-nod positions; underestimating the 
corrections will result in artificially steep spectral index 
estimates. In the top panel of FigurelHwe also plot a one- 



to-one line, showing that these two estimates are gener- 
ally consistent with one another. The dispersion about 
an ordinary least squares fit line is less than «0.05dex 
(i.e., 12%). 

In the bottom panel of Figure SI we plot the ratio of 
these two star-formation rate estimates, illustrating the 
effect for different assumptions of a^"^ . We allow a^"^ 
to range between 0.6 — 1.05 (left to right) in increments 
of 0.05. Except for the six discrepant sources identi- 
fied above, the two estimates appear generally consistent 
within errors when varying the assumption for the non- 
thermal radio spectral index. Given the extremely good 
agreement between these two estimates, which may not 
be terribly surprising given that the 33 GHz flux densities 
are generally dominated by thermal emission, star for- 
mation rates estimated using Equation [TUl arc assumed 
as our reference value for comparisons throughout the 
paper. By using Equation |TU1 we are able to calculate 
star formation rates for the entire sample (i.e., lower fre- 
quency radio data are not necessarily required). Before 
moving on, it is worth investigating how the variations 
in the physical size of the GBT beam may affect the star 
formation rate estimates given that EquationjTUldoes not 
take into account morphological differences between the 
(compact) free-free and (diffuse) non-thermal emission in 
galaxies. 

4.1.1. Distance Effects 

The distances to the galaxies being investigated vary 
by a factor of '^10. Consequently, the physical areas our 
observations average over span a factor of ^^100. By aver- 
aging over larger pieces of each galaxy, more diffuse emis- 
sion is introduced into each beam relative to additional 
emission from compact Hii regions and, for the 33 GHz 
observations, likely introduces more non-thermal emis- 
sion. To investigate the consequence of this effect, we 
plot the 33 GHz thermal fraction estimates as a function 
of the projected diameter of the region being measured 
in Figure |3l We find large (i.e., w 93l]^Q %) thermal frac- 
tions for all sources having sizes <0.5kpc, and a larger 
mix of thermal fractions for larger diameters. Excluding 
the 6 sources having thermal fractions <40%, as they 
may be artificially low (see above), there does appear to 
be a weak trend of decreasing thermal fraction with in- 
creasing projected area. However, it is hard to reliably 
quantify any trend given that there are so few sources 
with lower frequency radio data at the large distance end 
of the sample. 

To look at the effect of distance in another way, we plot 
the 33 GHz star formation rates, normalized by the pro- 
jected area of the beam, versus distance in Figure |5l for 
all detected, non-AGN sources. We calculate star forma- 
tion rates using Equation[TOl with a^'^ — 0.85, since this 
expression takes both thermal and non-thermal emission 
components into account. If the star formation rate cal- 
ibration is being affected by the inclusion of a substan- 
tial amount of non-thermal emission for the more distant 
galaxies, we would expect a trend of increasing star for- 
mation rate per unit area as the beam size increases. 
Perhaps there is a weak trend; the median star forma- 
tion rate surface density at distances below and above 

15Mpc is sa0.03 and «0.08 yr~-^ kpc^, respectively. 
However, the lack of a meaningful number of sources, es- 
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Figure 8. Top,- The 24 /im scaling factor for computing Ha + 
24 /im star formation rates plotted against the 33 GHz star for- 
mation rates for all non-AGN detected regions. Bottom: The IR 
scaling factor for computing FUV+IR star formation rates plotted 
against the 33 GHz star formation rates for all non-AGN detected 
regions. 

pecially extranuclear regions, at d > 15Mpc, makes it 
difScult to quantify any trend. Thus, is it appears that 
larger beams do not significantly affect our 33 GHz star 
formation rate estimates. 

4.2. Comparison with 84 fJ-m Warm Dust Emission 

Unlike the far-infrared continuum emission, for which 
a non-negligible fraction m ay be powered by an 
older stellar population (e.g. , iSauvage &: ThuanI Il992l : 
IWalterbos fc Greenawal^llQQa iBaldo et al.ll2010D . warm 
dust emission appears to be more tightly correlated to 
the curre nt star formation activity in the disks o f galax- 
ies (e.g.. iHelou eral|[200l fMurphv et al.ll2006aD . Con- 
sequently, a number of calibrations relating the warm 
dust emission at 24 /zm to star formation rates have 
been introduced in the lite r ature (e.g., IWu ct al. 2005; 
Perez- Gonzalez et al. 2006; Alonso-Hcrrcro ct al. 2006; 
Calzetti et al...2007; .Relafio et al, 2007: Zhu et al.. 2008; 



Rieke et al.ll2009() ~a detailed c omparison among eac h of 



these relations can be found in iCalzetti et al.l (|2010[ ). 

In Figure |6] we plot 24 fim flux densities against cor- 
responding (corrected) 33 GHz flux densities for the en- 
tire sample of 103 galaxy nuclei and extranuclear regions. 
Apart from the galaxy nuclei that are known to harbor 
AGN, a number of which are radio loud at 33 GHz rela- 
tive to the observed 24 fim emission, there is a fairly tight 
correlation between the 24 fim and 33 GHz flux densi- 
ties. The only other clearly outlying source is NGC 1377, 



which has a large 24 /im flux and remains undetected at 
33 GHz. This galaxy is thought to be experiencing a 
nascent starburst, for which it is within a few Myr of 
the onset of an intense star formation epi sode after be- 
ing quiescent for at least ~100Myr (e.g.. IRoussel et ahl 
12003^ .2006 ). This starburst, while heating the dust, has 
yet to produce detectable free-free emission, nor observ- 
able signatures of CRs. Over-plotted in Figure |6] are 
both a one-to-one line, scaled by the median 24 /im to 
33 GHz flux density ratio < /,y(24/im)/S'33GHz >= 142, 
and an ordinary least squares flt to all non-AGN sources 
detected at 33 GHz given by. 



/.(24A.m) 
mJy 



= 140 



'33 GHz 

mJy 



0.94 



(12) 



which are consistent within errors. The scatter about the 
ordinary least squares fit is ~0.26 dex. 

After converting the 24 fim flux densities into lumi- 
nosities, we relate these to corresponding 33 GHz star 
formation rate estimates to come up with the following 
relation among all non-AGN detected sources. 



SFR. 



24 /im 



Mgyr- 



-1 



3.1 X 10 



-38 



vL^{2A /im) 



erg s~ 



0.88 



(13) 



The fltting coefficients are consistent wit h others in the 
literature (for a detailed compilation see ICalzetti et al.l 
l2Q10f) . In Figure [7] we plot the ratio of the 24 /im to 
33 GHz star formation rate against the 33 GHz star for- 
mation rates finding that the scatter about this relation 
is Ri0.25dex. The square identifies the l ocation of the 
anomalous dust d etection in NGC 6946 (iMurphv et al.l 
[2010; Scai fe et al.l [20101. which may also explain why 
other sources lie significantly below the unity line. 

We might expect different relations between warm 
24 /im dust emission per unit star formation rate among 
the extranuclear regions and nuclei due to additional old 
stellar population heating of grains by galaxy bulges. To 
investigate this, we can separate our sample, and mea- 
sure the median ratio of the 24 /im spectral luminosity to 
33 GHz star formation rate for the nuclei and extranu- 
clear regions. The average ratios, along with the associ- 
ated dispersions, are given in Table[5l Within the scatter, 
which is just under a factor of ^2, the median ratios are 
consistent with one another, as well as with the median 
ratio for all non-AGN detected sources. Comparing the 
ratio of the total IR luminosity to 33 GHz star formation 
rate for nuclei, extranuclear regions, and all non-AGN de- 
tected sources, we again find that the median values are 
consistent. Accordingly, we can use these average ratios 
to write empirical relations between the 24 /tm spectral 
luminosity and total IR luminosity with the star forma- 
tion rate that should be reliable to within a factor of ^2 
such that, 



( = 2.45 X 10-« 



t^Ly(24 /im) 



erg s 



and 



/ SFR 



Avg. 
IR 



3.15 X 10 



-44 



L 



IR 



(14) 



(15) 



\ Mq yr~i / V*^''S ^ 

We note that the the coefficient equating the to 
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tal IR lumino s ity w ith star formation rate derived in 
iMurphv et all (|2011[ ) is «23% larger than this empiri- 
cally measured coefficient, and thus consistent given the 
nearly factor of ~2 scatter. 

4.3. Comparison with Ha and 24 fJ-ni Emission 

Given that not all of the UV/optical photons will be 
absorbed and re-radiated by dust, a series of new em- 
pirical calibrations based on the linear combination of 
observed 24 /itm (obscured star formation) and Ha (unob- 
scure d star for mation) lum inosities have been developed 
(e.g., i Calzetti et al. 2007; 'Ken nicutt et al.l 120071 120091: 
iZhu et al.,,2008 ). These empirical star formation recipes 
usually take the form of 



SFR„ 



M, 



= 5.37 X 10 



-42 



Lua +64i^^i/(24^m) 



erg s~^ 

(16) 

where the coefficient ^24, scaling the warm 24 /im dust 
emission, has been found to vary depending on the phys- 
ical scale being investigated. For in stance, on the scale o f 
individual Hii regions, ^24 « 0.031 (ICalzetti et al.ll2007n. 
wher e as for whole galaxies ^24 ~ 0.020 ()Kennicutt et al.l 
l2009f) . A luminosity dependence has also been found, 
with £24 increasing w ith the 24 fim luminosity of a galaxy 
(ICalzetti et al.ll2010() . 

Combining the 33 GHz star formation rates with the 
available Ha and 24 /im photometry, we can estimate 
£24 for the star-forming (i.e., non-AGN) regions in our 
sample. These values are plotted in the top panel of 
Figure |8] against the 33 GHz star formation rates, show- 
ing a large scatter. The median and dispersion are 
C24 = 0.018 ±0.004 and a^^^ = 0.031, respectively. Thus, 
on ^kpc scales, which likely average over many non- 
coeval Hii regions, the scaling coefficient appears to be 
similar to that for entire galaxies. It is also worth noting 
that the scatter is quite large, which most likely arises 
from both the large range in physical sizes being covered 
by the GBT beam among the entire sample, as well as the 
difficulty in matching the GBT and imaging photometry. 

4.4. Comparison with Total IR and UV Emission 

Similarly, by combining the star formation rate esti- 
mates from the total IR and (observed) UV emission, 
one can account for the obscured and unobscured emis- 
sion contributing to the total (bolometric) star formation 
rate, which we define as: 



(17) 

This diagnostic is often used to characterize star 
formation rates from ga laxies in b oth l ow- (e.g., 
Ilglesias-Paramo et al.l 120 06; Buat et al. 20071120111) and 



SFRtot = SFRpuv + SFRiR. 



high-z (e.g., [I glesias-Paramo et al. 2004 



. . lElbaz et al.l 

I2007t IDaddi et" al. 2007; Rcddy ct al. 201^ studies. 

Using the above calibration for the FUV star formation 
rate (i.e., Equation[5]), this relation can be expresses as a 
linear combination of the UV and total IR emission such 
that, 

-^^^'l = 4.42 X 10-44 ( ^FUV+6Ri:-iRA 
Moyr-iy V ergs-i ' 

where the coefficient £ir, scaling the total IR emission, 
has been found be ~ 0.46 for normal star-forming galax- 
ies (|Hao et al.ll2011[ ) and as large as ~0.6 for starburst 
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Figure 9. The ratio of tiie 33 GHz to total infrared (IR; 
8 — 1000 fim) fluxes plotted as a function of radio spectral index 
measured between 1.7 and 33 GHz. Nuclei and extranuclear regions 
are shown by circles and crosses, respectively, and AGN are iden- 
tified by filled circles. The solid line is the ordinary least squares 
fit to tfie data, while the dotted line illustrates the expected trend 
for a fixed ratio of t'Si.y ghz/^ir- We show the names for sources 
having both flux ratios i^Saa qhz/^ir > 4 x 10~^ and spectral 
indices of'^'Joll^ < 0.5 as these sources lie near the position of 
NGC6946 'Enuc.4 IIMurphv et al.|[20Tol : IHcaife et ahllSMoTl in this 
plot, and may be good candidates for containing anomalous dust 
emission. 



galaxies fe.g.. IMeurer et al.lll999l: ICalzettill2001[ ). These 
empirical values a re significantly smalle r than the value 
of 0.88 reported in IMurphv et al.l (|2011[ ). which assumed 
that the entire Balmer continuum was absorbed and rera- 
diated in the infrared for their derived IR star formation 
rate calibration of individual star-forming complexes. 

Using our 33 GHz star formation rates, we estimate £ir, 
for the star-forming regions in our sample and plot them 
in the bottom panel of Figure |S| against the 33 GHz star 
formation rates. The median and dispersion are S,ir — 
0.50 ± 0.05 and cr^jp — 0.47, respectively. The scatter 
is even larger than that for £24, which may arise from 
the fact that a calibration such as this depends on the 
range of ages among individual star formation sites in 
each star-forming complex covered by the GBT 33 GHz 
beam. Thus, we expect increasing uncertainties in such a 
calibration when applied to a single star-forming region 
depending on how well the IR and UV are measuring 
emission from recent star formation, not to mention the 
fact that the IR and UV are sensitive to star formation on 
different timescales (i.e., ~10Myr versus '^10 — 100 Myr, 
respectively) . 

5. DISCUSSION 

We have used new GBT 33 GHz radio continuum data 
to measure free-free emission toward 103 nearby galaxy 
nuclei and extranuclear star-forming regions. We find 
that the 33 GHz emission is typically dominated by free- 
free emission, making it a sensitive measure of the cur- 
rent star formation rate. Here we discuss the results 
of the comparisons between the 33 GHz star formation 
rates and other dust-sensitive estimators. We addition- 
ally discuss some instances where the 33 GHz emission 
may contain a contribution from anomalous dust emis- 
sion and how having high frequency radio data can help 
infer physical properties of the star-forming complexes. 
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Table 5 

Average Infrared to 33 GHz Star Formation Rate Ratios 



Sources < logpFRaa GHz/i^-J-i'(24 /xm)] > a < log(SFR33 ghz/-J-ir) > cr 
(dex) (dex) (dex) (dex) 



Nuc. -42.61 0.27 -43.55 0.28 

Enuc. -42.61 0.26 -43.49 0.29 

All -42.61 0.26 -43.50 0.29 



Note. — Median ratios arc given. Only non-AGN detected sources included. 
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Figure 10. The ratio of the 33 GHz to total IR fluxes plotted as 
a function of projected area of the GBT 33 GHz beam. Nuclei and 
extranuclear regions are shown by circles and crosses, respectively, 
and AGN are identified by filled circles. There does not appear to 
be a strong trend between the flux ratio and increasing projected 
area. 

5.1. Dust- Dependent Star Formation Rates 

Using the 33 GHz star formation rates as a reference, 
we have investigated a number of dust-sensitive star for- 
mation rate indicators used in the hterature. In tj4.2l we 
find a relation between the 33 GHz star formation rates 
and 24 /im luminosity that is consistent with others in 
the literature. However, while similar, we do find sig- 
nificantly larger scatter. We believe that the increased 
scatter most likely arises due to a combination of several 
effects. First, properly matching the single-beam GBT 
photometry to that from the ancillary images of resolved 
sources is difficult, and the associated uncertainties will 
increase the dispersion of any trends. Additionally, the 
area covered by GBT beam covers a range of physical 
sizes in the sample, which will also inject scatter into 
any observed correlations. And, of course, the correc- 
tions applied in cases of potential over subtraction of sky 
due to reference nods landing on bright regions of the 
targeted galaxy add another source of scatter. 

In looking at the relation between the 33 GHz star for- 
mation rates and 24 /xm spectral luminosities for galaxy 
nuclei and extranuclear regions independently, we find 
them to be consistent. This is also true for the rela- 
tion between 33 GHz star formation rates and total IR 
luminosities. Thus, on average, there does not appear to 
be excess infrared emission per unit star formation for 
galaxy nuclei relative to the extranuclear star-forming 
complexes on the physical scales being investigated. 

In i j4.3l and 14. 4[ we investigate two hybrid star for- 
mation diagnostics that attempt to account for the un- 



obscured and obscured star formation components (i.e., 
Ha + 24 /zm and UV-I-IR) . While we find scaling coeffi- 
cients for the 24 /^m and total IR luminosities that are 
similar to those reported in the literature, albeit with 
a large scatter most likely due to the reasons mentioned 
above, it is worth noting that the values reported are sim- 
ilar t o those found from measurements of e ntire galaxies 
(e.g.. lKennicutt et al.|[2009l : lHao et al.][20lll) . This agree- 
ment most likely arises because of the large area covered 
by the GBT beam at the distances of the sample galax- 
ies, being roughly '^Ikpc, on average. Thus, on these 
physical scales, and by not including a local background 
subtraction, it appears that the star formation diagnos- 
tics presented here are applicable for globally integrated 
measurements of galaxies. 

We summarize the final relations for each of these four 
dust-dependent star formation rate estimates, along with 
all theoretically motivated calibrations discussed in ij3.3l 
in Appendix Table |S1 

5.2. Variations on the Thermal Fraction and 
Non-Thermal Spectral Index 

In Figure [H we plot the ratio of the 33 GHz to total 
IR fluxes as a function of the radio spectral index mea- 
sured between 1.7 and 33 GHz. There is a clear trend 
of increasing flux ratio with flatter spectral index. This 
trend does not seem to arise due to differences in the 
projected physical area of the 33 GHz beam among the 
sources (see Figure fTU]) . In fact, this trend is opposite the 
expectation for distance effects; as the physical area sub- 
tended by the GBT beam increases, one would expect a 
larger contribution of diffuse non-thermal emission from 
a galaxy disk relative to additional free-free emission, 
thus resulting in a steeper spectral index and an increase 
in 33 GHz to IR flux ratios. 

The ordinary least squares fit to the non-AGN sources 
(solid line) appears generally consistent with, albeit 
slightly flatter than, the expected trend assuming a fixed 
ratio of t^5'i.7GHz/^iR; the dotted line is given by scaling 
the median ratio i'S'i.7ghz/-Fir by (33/1.7)-^"". By only 
having radio data at two well-spaced frequencies, it is not 
possible to distinguish whether the spectral fiattening is 
being driven by an increase in the thermal fraction or 
fiattening of the non-thermal radio spectral indices. If 
we are to assume that this trend is in fact dominated by 
an increase in the thermal fraction, we can use the scat- 
ter about the ordinary least squares fit to the data to set 
the maximum dispersion for the non-thermal spectral in- 
dices among the non-AGN sources, which is found to be 
CTcNT < 0.13. Interestingly, this is similar to the disper- 
sion in global non- thermal spectral indices measured by 
'Nikla s et all (ITOOTI CTc^NT = 0.13). 

Again, assuming that the trend is in fact driven by a 
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change in thermal fraction among the sample, suggests 
that the non-thermal emission component may be quite 
similar among each star-forming region, with a nearly 
uniform CR electron injection spectrum and similar spec- 
tral steepening from associated energy losses (e.g., syn- 
chrotron, inverse Compton, bremsstrahlung, and ioniza- 
tion processes). More specifically, there must only be 
moderate variability in the ratio of synchrotron to total 
CR electron cooling processes among each of the star- 
forming complexes being investigated. 

5.3. Anomalous Dust Candidates 

Cosmic Microwave Background experiments were the 
first to discover the presence of an "anomalous" dust- 
correlated emission c omponent at frequ encies between 
~ 10 - 90 GHz (e.g., iLeitch et all [19971). In the origi - 
nal GBT pilot study of NGC 6946. IMurphv et all (pOlO ) 
identified an extranuclear star-forming complex that 
had a significant amount of excess emission at 33 GHz 
relative to what was expected by extrapolating from 
multifrequency measurements at lower frequencies (i.e., 
< 10 GHz). While anomalous dust emission has been ob- 
served in excess of synchrotr on and free-free compo nents 
in Galactic Hil regions fe.g.. IDickinson et al.ll2009() . this 
observation of excess 33 GHz emission towards a star- 
forming complex in NGC 6946 is identified as the first de- 
tection of anomalous dust emission outside of the Milky 
Way. Now, using the full sample data, we can look for 
other potential extragalactic anomalous dust emitting 
candidates. 

As stated above. Figure [S] shows a clear trend in 
which the ratio of 33 GHz to total IR fluxes increases 
as the radio spectral index measured between 1.7 and 
33 GHz flattens. In Figure [9] we also show the location 
of NGC 6946 E nuc. 4, identified as co ntaining anomalous 
dust emission ([Murphv et al.l I2010t Scaifc c Fall I2O1O0 
that accounts for «50% of the observed 33 GHz flux den- 
sity. However, this source does not appear to be signifi- 
cantly discrepant from the expected trend in this plot. It 
therefore appears that identifying anomalous dust candi- 
dates with such coarse radio spectral resolution may not 
be possible, and a much finer sampling (i.e., better than 
a factor of 2 in frequency steps) is necessary. 

Besides NGC 6946 Enuc. 4, there are a number of other 
sources that both show flat radio spectral indices with 
large 33 GHz to infrared flux ratios; sources having flux 
ratios i^S'sa ghz/^ir > 4 x 10~^ and observed radio spec- 
tral indices ap-^Q^^ < 0-5 have been highlighted as po- 
tential anomalous dust candidates in Figure [9] However, 
a detailed investigation of these regions to confirm or 
deny the presence of anomalous dust emission is outside 
the scope of this paper, and really requires additional 
radio data at finely spaced frequencies much closer to 
33 GHz (e.g., between ^^15 and 90 GHz) to confirm a 
peaked spectrum. 

5.4. The Nascent Starburst NGC 1377 

The archetypal nascent starburst NGC 1377, which is 
thought to be within a few Myr of the onset of an in- 
tense star formation episode after being quiescent for at 
least ~ 100 Myr, remains undetected in the radio even at 
the depth of our 33 GHz observations. Given its mea- 
sured flux density at 24 /zm and the 3 a upper limit of 



0.44 mJy at 33 GHz, we would have exp ected a ~70g 
detect ion. As has already been shown by iRoussel et al.l 
((200l . this non-detection in the radio most likely does 
not arise from the source being optically thick. Using 
the 33 GHz upper limit, we can set even more stringent 
constraints on this. Assuming an electron temperature 
of ~10^K and tssghz ~ 1 sets the expected brightness 
temperature and emission measure of the source to be 
~6000K and ~ 5 x lO^pccm"^, respectively. The cor- 
responding size of the source would therefore have to 
be <lpc with an electron density of > 7 x lO^cm"^. 
These values are significantly more extreme than the size 
and elect ron density lim i ts of < 28pc and >2700cm~^ in- 
ferred bv IRoussel et al.l ()2003[ ). 

Given the apparent absence of AGN activity in 
NGC 1377, any free- free emission associated with ongo- 
ing star formation powering the bright dust continuum 
appears significantly s uppressed. If the nasc ent starburst 
scenario suggested by IRoussel et al.l ()2003[ ) is true, the 
lack of free-free emission detected at 33 GHz could very 
well arise from having extremely young, deeply embed- 
ding star formation in which the dust is absorbing ^95% 
of the ionizing photons. The absorption of such a high 
fraction of ionizing photons by dust would then lead to 
star formation rates that are overestimated by a factor 
of ^2 by using calibrations relying on dust emission. 

6. CONCLUSIONS 

In the paper we have presented GBT 33 GHz photom- 
etry for 103 galaxy nuclei and extranuclear star-forming 
regions included as part of the Star Formation in Radio 
Survey. The sample galaxies are included in the Spitzer- 
SINGS and iJersc/ie/-KlNGFISH legacy programs, and 
therefore have large amounts of ancillary data for future 
followup, including mid- and far-infrared spectral map- 
pings. Here we summarize our main conclusions from 
this initial investigation. 

1. Among the non-AGN sources having lower fre- 
quency radio data, we find a median thermal frac- 
tion at 33 GHz of «76% with a dispersion of «24%. 
For all sources resolved on scales ^0.5 kpc, the ther- 
mal fraction is even larger, being >90%, on aver- 
age, however, there is very little lower frequency ra- 
dio data for galaxies at the large distance end of the 
sample. This suggests that the rest-frame 33 GHz 
emission provides a sensitive measure of the ioniz- 
ing photon rate from young star-forming regions, 
thus making it a robust star formation rate indica- 
tor, especially at high spatial resolution. 

2. We find an increase in the 33 GHz to total IR flux 
ratios as the radio spectral index flattens, which 
does not appear to correlate with the projected 
area of the GBT beam, but is instead consistent 
with the expectation from an increase in the ther- 
mal fraction among the sources. Consequently, 
the ratio of non-thermal to total IR emission ap- 
pears relatively constant among star-forming re- 
gions, suggesting only moderate variability in the 
CR electron injection spectrum and ratio of syn- 
chrotron to total cooling processes. Given the 
scatter in this trend, the maximum dispersion in 
the non-thermal radio spectral indices among these 
non-AGN sources is (t„nt < 0.13. 
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3. Using the 33 GHz star formation rates as a refer- 
ence, we derive scaling coefficients for the following 
recipes that attempt to sum the obscured and un- 
obscurcd star formation components: Ha + 24 /zm 
and UV+IR. Although we have targeted galaxy nu- 
clei and extranuclear star- forming complexes, with 
a median resolution of '^Ikpc, the hybrid scaling 
coefficients are consistent with those in the litera- 
ture derived for globally integrated measurements 
of galaxies. 

4. Identifying anomalous dust candidates based on a 
coarse sampling of the radio spectrum, even with 
a large lever- arm spanning 1.7, and 33 GHz, is in- 
conclusive. A much finer (i.e., better than a factor 
of 2) sampling in frequency space, spanning ~15 
to 90 GHz to actually measure the peak, appears 
necessary for conclusive detections. 

5. Given the depth of our 33 GHz photometry, we are 

able to put greater doubt on the possibility that 
NGC 1377 is radio faint due to being optically-thick 
at radio frequencies; the source size would need to 
be <1 pc with an electron density of > 7x 10^ cm~^. 
Assuming that the source is at the onset of a star- 
burst, the lack of detectable free- free emission at 
33 GHz would require dust to absorb >95% of the 
ionizing photons, leading infrared-derived star for- 
mation rates to be overestimated by as much as a 
factor of 2. 

We are currently in the process of collecting lf4-sized 

33 GHz maps using the VLA in the D-configuration for 
each of the sources targeted with the GET, along with a 
few others at lower declination. These maps will have an 
angular resolution of ~2", allowing us to sample physical 
scales ranging between 30 — 300 pc given the distances of 
the galaxy sample. At this resolution, we will be able 
to probe the scales of individual giant Hll regions, where 
we expect even larger thermal fractions for more accu- 
rate estimates of star formation rates. These maps will 
additionally enable us to isolate the location of potential 
extragalactic anomalous dust emitting sources with sig- 
nificantly greater accuracy. An initial round of data has 
been taken between October 2011 and January 2012. 
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APPENDIX 

Here we present the more detailed results from the GBT observations, as well as summarize all star formation rate 
equations either utilized or empirically derived in the present analysis. For the GBT observations, we give the flux 
densities and errors reported from the four individual CCB ports, and provide details on how much time was spent 
on each source in Tables IH] and [71 We again note that the flux densities of the individual channels reported here do 
not include any corrections. Upper limits for the individual channels are not given, but rather actual measured values 
and estimated errors are given. 

In Table |S] we provide a summary of the theoretical and empirically derived star formation rate calibrations from 
this paper. A thorough discussion of the theoretical derivations can be found in i )3.3l and in iMurphv et alj ( 20111) . 



while details on the empirically derived, dust-dependent star formation rate estimates are given in §4.21 and 14.31 We 
also provide the timescale for the emission of each star formation rate diagnostic to fade to 50 and 5% of its peak value 
after having had 100 Myr of continuous star formation come t o a halt. These timescales, which are sensitive to the 
exact star formation history, were estimated using Starburst99 (jLeitherer et al.|[T999i ) following the same assumptions 
that went into deriving the corresponding theoretical relations as described in ^3.3l and in the table caption. Additional 
considerations for the radiating times of CR electrons were included when calculating the fading time for the non- 
thermal radio continuum emission, the details of which can be found in iMurphvl (|2009f f. 
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Nuclear Photometry for Individual CCB Ports 
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U.D ± U.uy 


33, 


.31 


21/30 


14.0 


In VjrU 4oDy 


o Q _l_ n oc 


n _i_ n 01 
2.U it U.2i 


1.8 ± 0.12 


1.5 ± 0.20 


33, 


.58 


A I A 

4/ 4 


T 
2.7 


MO O /I cTn 
In vjU 4o7y 


ot: 1 _l_ r\ n A 


OA '7 _L n on 
24.7 ± U.2U 


OQ Q _L n 1 

2o.o ± u.12 


00 a _l_ [\ 00 
22. D ih U.22 


32, 


.98 


A 1 A 

4/ 4 


T 
2.7 


MO o /I Kn/t 
IN VjrU 4ay4 


(O.O it U.Uo 


'TA c _i_ n n« 
( 4.0 it U.UO 


( i.o it u.u ( 


DO. 4 it U.i4 


32 


.95 


1 /t /I /I 
14/14 


n Q 

y.o 


MOO ACf^n 


U.O ± U.U4 


n /I _i_ n no 
U.4 ± U.Uo 


n 1 n no 

0.3 ± U.Uo 


U.O ± U.U4 


33 


.21 


50/50 


33.3 


NGC4631 


6.4±0.14 


4.7ib0.15 


3.4±0.11 


3.0 ±0.20 


34 


.27 


4/ 4 


2.7 


NGC 4725 


0.3 ±0.05 


0.3 ±0.04 


0.3 ±0.06 


0.2 ±0.07 


32 


.01 


42/74 


28.0 


NGC4736 


3.7±0.15 


3.2ib0.15 


2.6±0.11 


2.5 ±0.21 


33 


.26 


4/ 4 


2.7 


NGC 4826 


7.7±0.18 


6.4ib0.16 


5.7±0.12 


5.1 ±0.20 


33 


.61 


4/ 4 


2.7 


NGC 5055 


3.2 ±0.47 


2.9ib0.37 


2.4 ±0.33 


2.0 ±0.76 


32 


.37 


4/ 4 


2.7 


NGC 5194 


7.7 ±0.45 


7.6 ±0.39 


5.6±0.31 


5.4 ±0.72 


32 


.91 


4/ 4 


2.7 


NGC 5457 


-0.0 ±0.44 


-2.2 ±0.39 


-1.6±0.46 


-0.9± 1.11 


31 


.05 


4/ 4 


2.7 


NGC 5474 


0.3±0.12 


0.1 ±0.08 


0.1 ±0.08 


0.2±0.11 


32 


.99 


53/57 


35.3 


NGC 5713 


8.5 ±0.43 


7.2 ±0.37 


6.1±0.14 


5.3 ±0.28 


33 


.93 


4/ 4 


2.7 


NGC 5866 


6.7±0.23 


6.3±0.17 


6.4 ±0.08 


6.3 ±0.16 


32 


.77 


12/12 


8.0 


NGC 6946 


19.9 ±0.16 


17.4±0.23 


15.3 ±0.18 


13.9 ±0.25 


33, 


.55 


4/ 4 


2.7 


NGC 7331 


2.2 ±0.22 


1.6 ±0.30 


1.3±0.19 


1.3 ±0.22 


34.48 


3/4 


2.0 
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Table 7 

Extranuclear Photometry for Individual CCB ports 



ID 




522.75 GHz 


'531.25 GHz 


'S'34.75 GHz 


'S'38.25 GHz 




-^nods 










(mJy) 


(mJy) 


(mJy) 


(mJy) 


(GHz) 


used/taken 


(min) 


NGC 628 Enuc. 1 


U.o ± u.oy 


n ct J_ n on 

U.b ± u.zy 


n a J_ n qt 
U.b ± U.o7 


n Q _L n Q 
U.o ± U.oo 


32.48 


7/ 7 


4.7 


NGC 628 Enuc. 2 


U.4 it U. 11 


y).o ± U.Uo 


n /I _L n nc; 
U.4 it U.UO 


n Q J- n nn 

U.o it U.uy 


QQ f;i 
vSo.Ol 


ACX / An 
4U/ 4U 


Ofi v 
2d. / 


NGC 628 Enuc. 3 


U.b ± U.12 


n cr 1 n no 

U.O ± U.Uo 


0.5 ± 0.05 


n ,1 1 n nn 

U.4 ± U.uy 


33.49 


37/40 


24.7 


NGC 628 Enuc. 4 


0.1 di 0.15 


n 1 _L n HQ 
U.l ± U.Uo 


0.1 ± 0.10 


n _L n 1 
U.2 it U. 12 


32.95 


37/60 


24.7 


NGC 2403 


Enuc. 


1 


l.y ± U.OD 


l.y ih U . Z ( 


1.6 ± 0. 14 


1.5 ± 0.30 


33.50 


4/ 4 


T 
2. / 


NGC 2403 


Enuc. 


2 


1.7 lb 0.33 


1.2 zt 0.24 


1.5 ± 0.14 


1.1 ± 0.29 


33.31 


4/ 4 


2.7 


NGC 2403 


Enuc. 


3 


4.2 ± U.oo 


3.6 =b 0.23 


Q n _L n 1 Q 
o.U ± U.lo 


Q _L n 00 
2.0 ± U.2o 


33.73 


4/ 4 


T 

2.7 


NGC 2403 


Enuc. 


4 


1.7 lb 0.12 


1.0 ± U.Uo 


1.3 ± 0.05 


1.2 ± 0.11 


33.37 


24/24 


16.0 


NGC 2403 


Enuc. 


5 


l.O ± U.oZ 


I. ( ± U. /U 


1 Q _L n /I 
1.0 ± U.24 


n n _L n /1 
U.y ± U.42 


34.62 


2/ 2 


1.3 


NGC 2403 


Enuc. 


6 


r\ n _L n 1 c; 

u.y ± U.lO 


U. i ± U.Uo 


n _L n nn 

U.o ± U.uy 


0.7 ± 0.11 


33.08 


59/59 


39.3 


NGC 2976 


Enuc. 


1 


4.6 lb 0.31 


4.U it U.z4 


3.4 lb 0.13 


Q _L n 0(^ 

0.2 ± U.2b 


33.46 


4/ 4 


7 
2. / 


NGC 2976 


Enuc. 


2 


Z.yj ± U.ol 


2.U ± U.zo 


2.1 lb 0.14 


1 T _L n 0'7 

1.7 ± U.27 


33.21 


A / A 

4/ 4 


T 
2.7 


NGC 3521 


Enuc. 


1 


U.o ± U.Uo 


U.2 ± U.U4 


n Q _i_ n n/1 
U.o ± U.U4 


n _L n nc 
U.2 ± U.Uo 


32.39 


65/66 


43.3 


NGC 3521 


Enuc. 


2 


1 o 1 n 1 Q 
l.z it U.lo 


1 Q _L n oi 
l.O it U.zl 


n n _i_ n m 

U.y it u.iu 


n c _L n 1 Q 
U.o it U.lo 


QQ QT 
00.0 ( 


4/ 4 


7 
2. r 


NGC 3521 


Enuc. 


3 


0.5 lb 0.05 


n /I _L n nc 
U.4 ± U.Uo 


n Q _i_ n nc 
U.O ± U.Uo 


n Q _L n nc 
U.O ± U.Uo 


32.77 


59/60 


39.3 


NGC 3627 


Enuc. 


1 


O.U ± u.iy 


2.5 lb 0.18 


2.3 lb 0.11 


_L n on 
2.2 it U.2U 


33.17 


A / A 

4/ 4 


'7 

2.7 


NGC 3627 


Enuc. 


2 


'7 A 1 ft 1 Q 

/ .4 it U. IcS 


O. / it U.lo 


c fi _i_ n 11 

o.y it u.ii 


c f; _L n in 

o.b it u. ly 


QQ Q1 

00. Ol 


4/ 4 


'7 
2. ( 


NGC 3627 


Enuc. 


3 


U.y ± U.lo 


n o _L n 1 n 

U.o ± u.iy 


0.5 lb 0.11 


n c _L n on 
U.O ± U.2U 


33.17 


4/ 4 


T 
2.7 


NGC 3938 


Enuc. 


1 


U.2 ± U.Uo 


n 1 _L r\ r\A 
U.l ± U.U4 


n 1 _L n no 
U.l ± U.Uo 


n 1 -I- n n IX 
U.l dz U.Uo 


33.01 


49/50 


32.7 


NGC 3938 


Enuc. 


2 


U.4 ± U.oo 


n A _L n /I /I 
U.4 it U.44 


n Q _L n c; Q 
U.O ± U.OO 


n Q _J_ n fiO 

u.o ± U.b2 


31.23 


44/ 44 


29.3 


NGC 4254 


Enuc. 


1 


1.4 dz 0.10 


1.2 zb 0.10 


1.0 lb 0.13 


n _L n 00 
U.o ± U.22 


31.92 


8/ 8 


5.3 


NGC 4254 


Enuc. 


2 


U.U ± U.lz 


n n -1- n n^ 
-U.U ± U.U/ 


n 1 _L n no 
-U.l ± U.Uo 


n n _L n 1 n 
U.U ± U.IU 


32.94 


50/50 


33.3 


NGC 4321 


Enuc. 


1 


0.1 lb 0.10 


n 1 _L n nn 

U.l ± u.uy 


n 1 _L n no 
U.l ± U.Uo 


-0.0 ±0.11 


32.52 


43/ 5 1 


28.7 


NGC 4321 


Enuc. 


2 


-0.4 lb 0.11 


n Q J_ n n^ 
-U.o ± U.U/ 


n _L n nv 
-U.2 ± U.U / 


n _L n nn 

-U.2 ± U.uy 


33.39 


49/49 


32.7 


NGC 4631 


Enuc. 


1 


1.1 lb 0.16 


1.0 lb 0.10 


1.0 lb 0.15 


n T 1 n 1 T 

U.7 ± U.17 


32.43 


25/50 


16.7 


NGC 4631 


Enuc. 


2 


U.y ± U.lz 


n o _L n nn 

U.o ± U.uy 


n T _L n nv 
U. / ± U.U / 


n c; _L n nn 

U.O ± U.uy 


33.62 


24/26 


16.0 


NGC 4736 


Enuc. 


1 


3.3 lb 0.15 


2.8 lb 0.15 


a _L n 1 
2.b it U.12 


2.1 ± 0.21 


33.15 


4/ 4 


T 

2. /^ 


NGC 5055 


Enuc. 


1 


U.4 ± U.21 


n 1 1 r\ o o 

U.l ± u.oo 


U.l ± U.2o 


no 1 n ,1 

U.2 ± U.24 


32.28 


17/17 


11.3 


NGC 5194 


Enuc. 


1 


U.o ± U.OO 


A Q _L n CO 

-4.0 ± U.OZ 


-0.0 ± U.OO 


n Q _L f\ Q A 

-U.O ± U.o4 


31.89 


Q / A 

3/ 4 


2.0 


NGC 5194 


Enuc. 


2 


o n j_ n ACi 

z.y it u.4y 


Q C -1- n At 
o.O it U.4o 


n n -1- n tA 
U.y it U.o4 


-U.l it U.o2 


QO Tl 
o2. / 1 


Q / Q 
0/ 


n 
2.U 


NGC 5194 


Enuc. 


3 


1.2 lb 0.46 


1.1 lb 0.47 


1.1 lb 0.35 


0.5 ± 0.75 


32.15 


4/ 4 


2.7 


NGC 5194 


Enuc. 


4 


1.7 ± U.4o 


O _!_ n AC 

2.7 ± U.4o 


1.0 ± U.o4 


C _L n TO 

2.0 ± U.7o 


32.22 


4/ 4 


7 
2.7 


NGC 5194 


Enuc. 


5 


1 a _L r\ Act 
l.y ± U.4d 


2.3 lb 0.51 


1 n _i_ n Qo 
1.9 ± U.o2 


1 _L n ci(\ 
1.2 ± U.D9 


32.39 


4/ 4 


T 
2.7 


NGC 5194 


Enuc. 


6 


1.0 lb 0.14 


n n _L n 1 o 
U.y ± U.12 


nolo nT 

U.O ± U.U7 


r> 1 r\ 1 

U.y ± u.12 


33.30 


23/23 


15.3 


NGC 5194 


Enuc. 


7 


U.o ± U.4o 


0.2 lb 0.51 


n ^; _L n qo 
U.b ± U.o2 


n T _L n £jn 

U.7 ± u.by 


32.41 


4/ 4 


7 
2.7 


NGC 5194 


Enuc. 


8 


1.1 lb 0.51 


1 c _L n ctr\ 
-1.0 ± U.oU 


1 _L n /I n 
2.1 ± U.4U 


1.8 ± 0.81 


32.05 


3/ 4 


2.0 


NGC 5194 


Enuc. 


9 


1 O _L r\ A A 

l.z it U.44 


n o _L n /I o 
U.o ± U.4o 


n (; _L n Q 
U.b ± U.OO 


n c; _L n vq 

U.O ± U. /o 


32.19 


4/ 4 


7 
2. / 


NGC 5194 


Enuc. 


10 


U.O ± U.4o 


1 O _L r\ A'7 

1.2 ± U.4 / 


n Q _L r\ A 
U.O ± U.o4 


n _L n 70 
-U.o ± U. /2 


32.23 


4/ 4 


2.7 


NGC 5194 


Enuc. 


11 


n _L n on 
-U.O ± U.zU 


n _L n o Q 
-U.O ± U.2o 


n Q _L n 1 f; 
-U.O ± U.lb 


n Q _L n /I 
-U.O ± U.24 


32.65 


14/14 


9.3 


NGC 5457 


Enuc. 


1 


n Q _L n 17 
U.o it U. 1 / 


n v _L n Qn 
U. / it U.oy 


n A _L n 00 
U.b it U.2o 


U.O it U.20 


Q1 OA 
o1.2D 


1 v / 1 7 
1 / / 1 / 


1 1 Q 
11.0 


NGC 5457 


Enuc. 


2 


0.3 lb 0.15 


n /I _L n 1 c 
U.4 ± U.lo 


n _L n 1 
U.2 ± U.lo 


n _L n 1 n 

U.2 ± u.iy 


31.93 


32/40 


21.3 


NGC 5457 


Enuc. 


3 


o.O it U.04 


o 1 1 n 4 c 
y.l it U.40 


n _L n c c; 
y.z it U.OO 


V /I _L 1 1/1 
/ .4 it 1. 14 


Q1 QA 
0l.04 


A / A 
4/ 4 


7 
2. / 


NGC 5457 


Enuc. 


4 


Z.4 ± u.oo 


Q Q _L n Of; 
O.O ± U.2d 


1 Q _L n Q v 

1.0 ± u.o / 


1.1 ±0.86 


31.48 


7/ 7 


A 7 

4.7 


NGC 5457 


Enuc. 


5 


z.o ± u.4/ 


o o _L n Qf; 
2.0 ± U.OO 


_L r\ AC\ 

2.0 ± u.4y 


2.8 ±1.12 


31.20 


4/ 4 


7 

2. / 


NGC 5457 


Enuc. 


6 


3.1 lb 0.52 


4.1 lb 0.38 


Q "7 _L C\ a. A 

0. / ± U.b4 


2.4 ± 1.37 


31.02 


3/ 3 


2.0 


NGC 5457 


Enuc. 


7 


5.0 lb 0.45 


5.1 lb 0.35 


4.8 lb 0.54 


4.5 ±1.21 


31.33 


4/ 4 


2.7 


NGC 5713 


Enuc. 


1 


1.8 lb 0.42 


1.8 lb 0.42 


l.OibO.16 


1.1 ±0.28 


34.28 


4/ 4 


2.7 


NGC 5713 


Enuc. 


2 


2.4 lb 0.42 


1.7 lb 0.38 


1.3ib0.16 


1.3 ±0.29 


34.08 


4/ 4 


2.7 


NGC 6946 


Enuc. 


1 


0.3 lb 0.09 


0.5 lb 0.08 


0.5ib0.10 


0.5 ±0.12 


31.70 


16/16 


10.7 


NGC 6946 


Enuc. 


2 


2.8 lb 0.17 


2.6 lb 0.22 


2.2 lb 0.21 


2.1 ±0.25 


32.64 


4/ 4 


2.7 


NGC 6946 


Enuc. 


3 


1.2 lb 0.10 


1.1 lb 0.09 


0.9ib0.12 


0.9 ±0.13 


32.52 


15/16 


10.0 


NGC 6946 


Enuc. 


4 


3.4ib0.10 


3.2 lb 0.10 


2.8 lb 0.12 


2.7±0.13 


33.07 


15/16 


10.0 


NGC 6946 


Enuc. 


5 


0.6 lb 0.07 


0.6 lb 0.07 


0.4 ±0.09 


0.5 ±0.08 


32.37 


37/37 


24.7 


NGC 6946 


Enuc. 


6 


3.6 lb 0.16 


3.1 ±0.22 


3.0 ±0.21 


2.2±0.26 


32.98 


4/ 4 


2.7 


NGC 6946 


Enuc. 


7 


3.6ib0.21 


3.5 ±0.24 


3.0 ±0.31 


3.1 ±0.29 


32.48 


3/ 4 


2.0 


NGC 6946 


Enuc. 


8 


2.1ib0.19 


1.9 lb 0.20 


1.7±0.26 


1.7 ±0.24 


32.41 


4/ 4 


2.7 


NGC 6946 


Enuc. 


9 


3.1 lb 0.20 


2.6 lb 0.20 


2.4 ±0.28 


2.2 ±0.23 


32.81 


4/ 4 


2.7 
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Table 8 

Summary of Star Formation Rate Conversions 



Observable Equation Fading Timcscale^ 



Theoretical Relations'^ 



Extinction corrected FUV flux 


SFR = 


4.42 X 10-'"'Lfuv 


~ 5, 100 Myr 


Extinction corrected NUV flux 


SFR = 


7.15 X 10-44LNUV 


~ 5, 150 Myr 


Extinction corrected Ha recombination line flux 


SFR = 


5.37 X 10-42 


~ 2, 5 Myr 


Free-free (radio) flux density 


SFR = 


4.6 X 10-28T-f«^o.i^iT 


~ 2, 5 Myr 


Non-thermal (radio) flux density 


SFR = 


6.6 X 10-29 i^gNT^NT 


~ 40, 100 Myr 


Total (radio) flux density 


SFR = 


10-2V(2.18TO-«^-0-^i + 15.1^-^^-)L, 


~ 2 - 40, 5 - 100 Myr 



Empirical Relations'^ 



24 /im flux density 


SFR = 2.45 X 10-"'='yLy(24/im) 


~ 5, 100 Myr 


Total infrared (IR; 8 - 1000 ^tm) flux 


SFR = 3.15 X 10-'''*LiR 


~ 5, 100 Myr 


24 fim flux density and Ha line flux 


SFR = 5.37 X 10-''2[Lhq + 0.018i/L^(24 ^tm)] 


~ 2 - 5,5 - 100 Myr 


Total IR and FUV fluxes 


SFR = 4.42 X 10-44 (Lpuv + O.SOLir) 


~ 5, 100 Myr 


^ The timcscalcs given arc those for the emission to drop to 50 and 5% of its peak value after having had 100 Myr of continuous 
star formation come to a halt. These values, which arc sensitive to the star formation history, were calculated using Starburst99 
ILeitherer et al.l|1999,) following the same assumptions that went into deriving the corresponding theoretical relations. The non- 
thermal radio continuum timescalc sums the fading timescalc of supcrnovae activity and appropriate fractions of the radiating lifetimes 



for 1.4 GHz emitting CR electrons in the ISM assuming a magnetic field strength of 5 fiG, equipartition between the magnetic and 
radiation field energy densities, and an ISM density of nisM — 0.1 cm~ (for details on the CR cooling timescales see lMurphvl[2009h . 
^ Each theoretical calibration was calculated using using Starburst99 dLeitherer et al.|199gD and assumes solar metallicity, continuous 
star formation, and a |Kroupa| j2001f) IMF, having a slope of —1.3 for stellar masses between 0.1 — 0.5 Mq and —2.3 for stellar masses 
ranging between 0.5 — 100 Mq. Star formation rates are in units of Mq yr~^ and observables are in cgs units. The UV calibrations 
are calculated using the GALEX FUV and NUV basspands, which are centered at 1528 and 2271 A, respecively. For radio-based star 
formation rates, i-'ghz and Te,4 are frequency and thermal electron temperature in units of GHz and 10'* K, respectively. Additionally, 

is the non-thermal radio spectral index, which is assumed to be ~ 0.85 in the present analysis. 
^ The empirical relations each have a scatter that is roughly a factor of <2. Consequently, the empirical calibration given here for 
the total IR luminosity is conistent with the theoretical relation given in IMurphv et~all 1 120111 . SFR = 3.88 x 10-''''Lir ), which, 
using starburst99 with the same assumptions above, assumed that the entire Balmer continuum was absorbed and reradiated in the 
infrared. See ^37311431 and 1431 for more details. 
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Figure 1. For each galaxy we show 2 side-by-side panels illustrating the GBT observations on 24 fim Spitzer images. The left p anel s 
indicate the location of each GBT 33 GHz pointing, while the right panels show the location of the GBT reference nod positions (see i|2.2l l. 
The diameter of each circle is 25", which corresponds to the typical FWHM of the GBT beam at 33 GHz. The 1 kpc linear scale bars 
assume the distances given in Table [T] 
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Figure 1. Continued 
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Figure 1. Continued 
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Figure 1. Continued 
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Figure 1. Continued 



